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Abstract 

The Averaged Null Energy Condition (ANEC) states that the integral along a complete null 
geodesic of the projection of the stress-energy tensor onto the tangent vector to the geodesic 
cannot be negative. ANEC can be used to rule out spacetimes with exotic phenomena, such as 
closed timelike curves, superluminal travel and wormholes. We prove that ANEC is obeyed by a 
minimally-coupled, free quantum scalar held on any achronal null geodesic (not two points can 
be connected with a timelike curve) surrounded by a tubular neighborhood whose curvature is 
produced by a classical source. To prove ANEC we use a null-projected quantum inequality, 
which provides constraints on how negative the weighted average of the renormalized stress- 
energy tensor of a quantum held can be. Starting with a general result of Fewster and Smith, 
we hrst derive a timelike projected quantum inequality for a minimally-coupled scalar held on 
hat spacetime with a background potential. Using that result we proceed to hnd the bound of 
a quantum inequality on a geodesic in a spacetime with small curvature, working to hrst order 
in the Ricci tensor and its derivatives. The last step is to derive a bound for the null-projected 
quantum inequality on a general timelike path. Finally we use that result to prove achronal 
ANEC in spacetimes with small curvature. 
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Chapter 1 
Introduction 


In the context of General Relativity, it is always possible to invent a spacetime with exotic 
features, such as wormholes, superluminal travel, or the construction of time machines, and 
then determine what stress-energy tensor is necessary to support the given spacetime according 
to Einstein’s Equations (units: h = c = 1) 

Gab = 9,T^Tab- ( 1 - 1 ) 


However, in quantum field theory, there are restrictions on Tab that could rule out exotic 
spacetimes. Two examples of these are quantum inequalities and energy conditions. 

Quantum inequalities (also called Quantum Energy Inequalities) are bounds on the weighted 
time averages of the stress-energy tensor. They were hrst introduced by Ford HZI to prevent 
the violation of the second law of thermodynamics. The general form of a quantum inequality 


IS 


dTf(t)T^(w(t))V‘V‘>-B 


( 1 . 2 ) 


where vj(t) is a timelike path parameterized by proper time t with tangent vector V, and / 
is a sampling function. The quantity R is a bound, depending on the function /, the path w 
and the quantum held of interest. 

Since then, they have been derived for a wide range of spacetimes, helds, and weighting 
functions. We concentrate here on quantum inequalities for minimally coupled scalar helds in 
curved spacetime. However, quantum inequalities for interacting helds have been derived in 
1-1-1 dimensions For systems with boundaries, there are diherence quantum inequalities 

[T8l [8], in which Tab in Eq. fll.2p is replaced by the diherence between Tab in some state of 
interest and Tab in a reference state. The bound B may also then depend on the reference 
state. However, such diherence inequalities cannot be used to rule out exotic spacetimes, at 
least in the case where the exotic matter that supports the spacetime comes from the vacuum 
state in the presence of the boundaries. 

Pointwise energy conditions bound the stress-energy tensor at individual spacetime points. 
One example is the Null Energy Condition (NEC) which requires that the null contracted 
stress-energy tensor cannot be negative. 


Tabtt>0, (1.3) 

for i°‘ a null vector. Classically, pointwise energy conditions seem reasonable, but in the 
quantum context they are violated. Quantum held theory allows arbitrary negative energy 
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densities at individual points, a well known example being the Casimir effect. Even in the 
simple case of a minimally coupled free scalar held, all known pointwise energy conditions fail 
and even local averages must admit negative expectation values [7]. 

On the other hand, averaged energy conditions bound the stress-energy tensor integrated 
along a complete geodesic; they are weaker and have been proven to hold in a variety of 
spacetimes. One example is the Averaged Null Energy Condition (ANEC) which bounds the 
null-projected stress-energy tensor integrated along a null geodesic 7 

TabttdX > 0 

To rule out exotic spacetimes such as those with wormholes and closed timelike curves, we 
would like to prove energy conditions that restrict the stress-energy tensor that might arise 
from quantum helds and show that the stress-energy necessary to support these spacetimes 
is impossible. We need a condition which is strong enough to rule out exotic cases, while 
simultaneously weak enough to be proven correct. 

The best possibility for such a condition seems to be the achronal ANEC [ 21 ], which requires 
that 7 of Eq. f|1.4p is a complete achronal null geodesic i.e., no two points of 7 can be connected 
by a timelike curve. That is to say, we require that the projection of the stress-energy tensor 
along a null geodesic integrate to a non-negative value, but only for geodesics that are achronal. 
As far as we know, there is no example of achronal ANEC violation in spacetimes satisfying 
Einstein’s equations with classical matter or free quantum helds as sourcesji] Achronal ANEC 
is sufficient to rule out many exotic spacetimes [ 2 T] . 

It has been proven that ANEC holds in Minkowski space and Ref. [n] showed that it also 
holds for geodesics traveling through empty, hat space, even if elsewhere in the spacetime there 
are boundaries or spacetime curvature, providing that these stay some minimum distance from 
the geodesic and do not ahect the causal structure of the spacetime near the geodesic. This 
proof made use of quantum inequalities for null contractions of the stress tensor averaged over 
timelike geodesics |14j . 

This work however, does not really addresses the possibility of exotic spacetimes. The 
quantum inequalities on which it depends apply only in hat spacetime, so they cannot be used 
to rule out spacetimes with exotic curvature. For that, we need limits on the stress-energy 
tensor in curved spacetimes, the work presented here. 

This thesis presents a complete proof of achronal ANEC for minimally coupled scalar helds 
in spacetimes with curvature in a classical background that obeys NEC, using a null projected 
quantum inequality. In the hrst chapter we present a general quantum inequality derived by 
Fewster and Smith [TB] that we use in later chapters to derive a bound. In the second chapter 
we derive the bound for hat spacetime with a background potential [25], a case similar to the 
curved spacetime one. In chapter three we present the timelike projected quantum inequality 
in curved spacetime [26] and discuss the importance of this result. In chapter four we use 
that result to derive a null projected quantum inequality, which we proceed to use to prove 
achronal ANEC [ 23 ] • Finally, in the Appendix we present a new class of coordinates, called 
multi-step Fermi coordinates and use them to write the connection and the metric in terms of 
the curvature [ 23 ], results we use throughout this work. 

We use the sign convention (—, —, —) in the classihcation of Misner, Thorne and Wheeler 
[28] . Indices a,b,c,... denote all spacetime coordinates while denote only spatial 

^Except for the case of non-minimally coupled quantum scalar fields. 



(1.4) 
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coordinates. We denote normal derivatives with comma and covariant derivatives with semi¬ 
colon. 



Chapter 2 

Absolute quantum energy inequality 


In this chapter we present a general quantum inequality derived by Fewster and Smith |15j . 
result which we are going to use throughout the thesis. First we consider a minimally-coupled 
scalar field with the usual classical stress-energy tensor, 

Tab = Va$Vb$ - • (2.1) 

where /i is the mass. Following Ref. |T5], we define the renormalized stress-energy tensor 

(TD = lim, rx"‘ mx)'P(x')) - H(x, x')) - . (2.2) 

. 7 ’— 


The quantities appear in Eq. fl2.2p will be defined below. is the point-split energy density 

operator, 

yspht = (g) ® Vrf/ -h • (2.3) 

The point-split energy density operator acts on the difference between the two-point function 
and the Hadamard series. 


H{x,x') = 


\_a+{x, x') 


71 + ^Vj{x,x')a^^{x,x') In 


3=0 


(^+{x,x') 

P 


+ Wj{x, x')a^(x, x') 

3=0 


(2.4) 


where a is the squared invariant length of the geodesic between x and x', negative for timelike 
separation. In fiat space 



cr(x, x') = —rjabix — x')“(x — x')’^. 

(2.5) 

By F(ct+), for some 

function F, we mean the distributional limit 



F((T+) = lim F((T,) , 

(2.6) 


e^0+ 

where 

a^ix, x') = a{x, x') + 2ie{t{x) — t{x')) -1- . 

(2.7) 


In some parts of the calculation it is possible to assume that the points have only timelike 
separation, so we define 

T = t-t' (2.8) 
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and write 


(2.9) 


F{<t+) = F(r_) = limF(T,), 

£—>■0 

where 

= T — ie . (2-10) 

We have introduced a length I so that the argument of the logarithm in Eq. fl2.4|) is dimension¬ 
less. The possibility of changing this scale creates an ambiguity in the dehnition of H, but this 
ambiguity for curved spacetime can be absorbed into the ambiguity involving local curvature 
terms discussed below. The ambiguity in the case of a field in flat spacetime with background 
potential is discussed in Ch. [3l For simplicity of notation, we will assume we are working in 
units where / = 1. 

The function A is the Van Vleck-Morette determinant bi-scalar, given by 


A(a:, x') 


det(Va <8) Vb'(j(a:, x')/2) 


( 2 . 11 ) 


The Hadamard coefficients are given by the Hadamard recursion relations, which are the 
solutions to 

{a +^‘^)H{x,x’) = 0. ( 2 . 12 ) 

The recursion relations for the minimally coupled field in a curved background are [T5] 

(□ -|- p )A ^ 2vQ^a(^’ T duo T vqOct = 0 , (2.13) 

(□ -|- /i^)nj + 2{j + l)vj+i^aO'’°‘ — 4:j{j + l)nj+i + {j + l)vj+iUa = 0 , (2-14) 

2^1,+ WiUa -\- 2ni_aCT’“ — 4ni -|- viUa = 0 , (2-15) 


(□ + + 2{j + - Aj{j + l)wj+i -F (j -h l)wj+ina 

+2vj^i^a<^’°‘ — 4(2j -|- l)nj+i + = 0 . (2-16) 

All the Vj, and the Wj for j > 1 are determined by the differential equations discussed above, 
but Wo is undetermined. Here we will follow Wald [iQ] and choose wq = 0. 

From Ref. [1^ we have the dehnition 

H(x, x') = ^ [H(x, x') + H(x', x) + iE(x, x')] , (2.17) 

where iE is the antisymmetric part of the two-point function. We can write Hj{x,x'), j = 
— 1, 0,1,..., to denote the term in H involving (with or without In (ct+)), and i7(j) to denote 
the sum of all terms up through Hj . We will split up E{x,x') into terms labeled Ej that are 
proportional to , dehne a “remainder term” 


and let 


j 

Rj = E — ^ ^ Ek , 

k=-l 


Hj{x, x') = - [Hj{x^ x') -|- Hj{x', x) + iEj{x, x'))] 
H(j){x, x') = ^ [H(j){x, x') + H^j){x', x) + iE{x, x'))] . 


(2.18) 
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(2.19a) 

(2.19b) 






The term Q in Eq. fl2.2p is the one introdnced by Wald to preserve the conservation of the 
stress-energy tensor. Wald [3S] calculated this term in the coincidence limit, 


Q = 


127r2 


Wi{x, x) 


( 2 . 20 ) 


The term Cab handles the ambiguities in the definition of the stress-energy tensor T in 
curved spacetime. We will adopt the axiomatic definition given by Wald [10], but there remains 
the ambiguity of adding local curvature terms with arbitrary coefficients. From Ref. [5] we 
find that these terms include 


^^'^Hab = 


^^^Hab = 


6 


6 




ab 


j yf^R^d^x = 2R, 
j ^R^^Rab = R-,ab - OR, 


ab - 2gabOR - + 2RRab 


2 

gabOR 


(2.21a) 


gabR R, 


■cd 




acbd • 


(2.21b) 


Thus in Eq. fl2.23l) we must include a term given by a linear combination of Eqs. fl2.21al) and 
fl2.21bp to first order in R, 

Cab = a^^^Hab + b^^^Hab ( 2 . 22 ) 

where a and b are undetermined constants]!] 

A spacetime is globally hyperbolic when it contains a Cauchy surface, a subset of spacetime 
which is intersected by every causal curve exactly once. Global hyperbolicity requires the 
existence of unique advanced and retarded Green functions. We define w{t), a timelike path 
contained in a globally hyperbolic convex @ normal neighborhood N and for this path we can 
state the quantum inequality of Ref. m- 


dtg{ty{TrV‘^V^)Mt) > 


d^ 


TT 


g » g(eXf"V‘V‘'H(^)) (-?.?) 


1 A 


dtg‘\t)i-Qgab + Cab)V‘^V’’. 


(2.23) 


The Fourier transform convention we use is 

f{k) or f^[k] = 


dxf{x)e 


ixk 


(2.24) 


In the inequality we take the transform with respect to both arguments and evaluate at ^ and 
The operator 6* denotes the pullback of the function to the geodesic, 

)((,«') = »((')), (2.25) 

^ There are also ambiguities corresponding to adding multiples of the metric and the Einstein tensor to 
the stress tensor. The first can be considered renormalization of the cosmological constant and the second 
renormalization of Newton’s constant. We will assume that these renormalization have been performed, and 
that the cosmological constant is considered part of the gravitational sector, so neither of these affects Tat- 

convex normal neighborhood N is one such that any point q € N can be connected to any other point 
p G N with a unique geodesic totally contained in N. For more detailed discussion of normal neighborhoods 
and their properties see Ref. [22] 
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and the subscript (5) means that we include only terms through j = 5 in the sums of Eq. fl2.4p . 
However, we will prove that terms of order j > 1 make no contribution to Eq. fl2.23p . 

We will use the general inequality of Eq. fl2.23p to provide a bound for the integral of the 
renormalized stress-energy tensor in three different cases. In Ch. [3] for the timelike projected 
T^i, in flat spacetime with a background potential, in Ch. 0] for the energy density in curved 
spacetime and in Ch. [5]for the null projected stress-energy tensor in spacetimes with curvature. 
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Chapter 3 

Quantum Inequality for a scalar field 
with a background potential 


As a first step toward deriving a bound for the quantum inequality in a spacetime with bounded 
curvature we first derive a quantum inequality in a flat spacetime with a background potential, 
i.e., a field with a mass depending on spacetime position. This is a simpler system that has 
many of the important features of quantum fields in curved spacetime. For a scalar field <F in 
a background potential, the Lagrangian is 


L = ^ - V{x)^‘^] , 

(3.1) 

the equation of motion is 


(□ -f l/(a;))$ = 0, 

(3.2) 

and the classical energy density is 


Tu = ^ + (V<l>)2 + U(a;)<F2] . 

(3.3) 

So the T^'"^ operator of Eq. (j2.3|) contracted with timelike vectors, 
background potential V becomes 

for a scalar field with 


^split _ 


^ dada' + 


V{x) + V{x') 


,a=0 


(3.4) 


where the potential is analogous to the mass of Eq. fl2.3p . Since we will take the limit where 
X and x' coincide the location of evaluation of V does not matter, but the form above will be 
convenient later. The renormalized stress-energy tensor in this case is 

= Im ((0(x)0(a;')) - H{x,x')) - Q{x '), (3.5) 

'T —V-r' 


As discussed in Chapter [2] there is an ambiguity in the above procedure. In order to take 
logarithms, we must divide a by the square of some length scale 1. Changing the scale to 
some other scale I' decreases H by 6H = 2{vo + Via -!-■■■) ln(/'//). This results in increasing 
Tab by \mix^x'{.dadh — {^/‘^)Vabd^dc)SH. Using the values for vq and Vi computed below, this 
becomes {l/12){y^ah — Thus we see that the definition of Tab must include 
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arbitrary multiple of {Vab — Vab^V). This ambiguity can also be understood as the possibility 
of including in the Lagrangian density a term of the form R{x)V{x), where R is the scalar 
curvature. Varying the metric to obtain Tat and then going to flat space yields the above term. 
The situation is very much analogous to the possible addition of terms of the form R^ and 
RabR°^^ in the case of a held in curved spacetime, which give rise to the local curvature terms 
in Eq. fl2.22p . 

Thus we rewrite Eq. (13.bh to include the ambiguous term, 

{nr(x')) = limy»“ mx)<l,(x')) - H(x,x')) - Q(x') + CU. . (3-6) 

x^x' 


where C is some constant. Whatever dehnition of one is trying to use, one can pick an 
arbitrary scale I and adjust C accordingly. 

So the quantum inequality of Eq. fl2.23p becomes 



dTg(tf{TS‘)(t,0)>-B, 


(3.7) 


where 


and 


B= r dtg^(t)(Q-CVu), 

J 0 J —oo 

F(t, (') = 9(*)s((')T““'“7?(5)((*. 0). (*'. 0)). 


(3.8) 

(3.9) 


F denotes the Fourier transform in both arguments according to Eq. fl2.24p . 

We work only in hrst order in V but don’t otherwise assume that it is small. We can 
express the maximum values of the background potential and its derivatives as 


\v\ < \VJ < 

lKal.|<C.x iKadidCL, 


where Vnax, Knax; Knax ^ax positive numbers, hnite but not necessarily small. 

First, we discuss the operator, then we compute the Hadamard series and we apply 
the operator. After that we perform the Fourier transform, leading to the hnal quantum 
inequality. 


3.1 General considerations 

We will now compute the quantum inequality bound B to hrst order in the potential V and 
its derivatives. In the next subsection, we will make some general remarks about terms in H 
coming from H, which are symmetrical under the exchange of x and x', and show that we need 
keep terms only through hrst order, not hfth order as in Eq. fl2.23p . Then we will simplify the 
operator dehned in Eq. fl3.4p . In the next section, we will compute, order by order, the 
terms and H. We will then take the Fourier transform to hnd the quantum inequality bound. 


3.1.1 Smooth, symmetrical contributions 


Dehne x = {x — x')/2, t = (t + t')/2 and t = t — t'. Let 


A(r) 



dt F 



(3.11) 


14 













so that = 74(—^). The presence of the g functions in Eq. fl3.9p makes F have compact 

support in t and t', so A has compact support in r. 

Suppose F contains some term / that is symmetrical in t and t'. Let a be the corresponding 
term in A according to Eq. flS.lip . Then a will be even in r, so a will be even also. If a G 
then a & and we can perform the integral of this term separately, giving an inverse Fourier 
transform, 

r f/(-?. ?) = r =«(o). (3.12) 

In particular, if 

lim/(f,f') =/W , (3.13) 

then 

/”-/(-?,0= r dtg(tf fit), (3.14) 

JQ J—oo 

and if f{t) = 0 there is no contribution. 

Terms arising from H appear symmetrically in FI. At orders j > 1 they have at least 
4 powers of r, so they vanish in the coincidence limit even when differentiated twice by the 
operators of Thus such terms make no contribution to Eq. fl4.4p . 


3.1.2 Simplification of 

We would like to write the operator of Eq. (13.41) in terms of separate derivatives on the 
center point x and the difference between the points. First we separate the derivatives in 
into time and space, 

3 

dada' = dtdf + Va; • Vx' ■ (3.15) 

a=0 

We can write the spatial derivative with respect tc0 x in terms of the derivatives at the 
endpoints, 

V'= V2+2V.-V.1 + V2,. (3.16) 

Then Eqs. (13.4p . (I3.15I) . (I3.16I) give 



1 

2 

1 

4 


3A' + 1 (V? - - V^,) + 1 (Vix) + V(x')) = 

[V? + Dx — d'l -\- Dx' — d'^t + 2dtdt' + V(x) + V (a^O] i 


(3.17) 


where and Dx' denote the D’Alembertian operator with respect to x and x'. Then using 


= :) “ 28 , 8 ,. + 8 ))] 


(3.18) 


we can write 




(□, + V{x)) H + (Dx^ + V{x')) H + VlH 


-dlH. 


(3.19) 


^When a derivative is with respect to t or t' (or x or x'), we mean to keep the other of these fixed, while 
when the derivative is with respect to t or r, we mean to keep the other of these fixed. When the derivative is 
with respect to x we mean to keep x — x' hxed. 
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Consider the first term. The function H{x,x') obeys the equation of motion in x, and so does 
E{x,x'). Thus all of H{x,x') is annihilated by + V{x), except for H{x',x), 

(□, + V{x)) H = ^{n, + V)H{x’, x) . (3.20) 

The quantities Vj{x, x') are symmetric, so the only sources of asymmetry in H are the functions 
Wj^ which are real but may not be symmetric, and the fact that the imaginary part of (t+{x, x') 
is antisymmetric. Thus 

H{x',x) = H{x,x'y + 'y^ywj{x', x) — Wj{x,x))ayx,x') , (3.21) 

j 

where * denotes complex conjugation. Since Da, is real, if we use Eq. fl3.2ip in Eq. fl3.20p . we 
have (Da, + V)H{x,x')* = 0, and we ignore Vwj, because it is second order in V, leaving 

(□a, + V{x)) H = -^Ox'^iwj{x',x) - Wj{x,x'))ayx,x'). (3.22) 

j 

By the same argument, 

(□a,/ + V{x')) H = 'Y^{wj{x,x') - Wj{x\x))ayx,x'). (3.23) 

j 

The hrst two terms in the brackets in Eq. fl3.19|) are the sum of Eqs. 03.221) . 03.23p . This sum 
is smooth, symmetric in x and x' , and vanishes in the coincidence limit. Thus according to 
the analysis of Sec. 13.1.11 it makes no contribution in the Fourier transform of Eq. 04.41) . and 
for our purposes we can take 

T^f^H 

3.2 Computation of H 

Examining Eq. 04.15P we see that is sufficient to compute H for purely temporal separation as 
a function of t, T, and x, the common spatial position of the points. The function H{t, t') is a 
series of terms with decreasing degree of singularity at coincidence: r“^, Inr, Inr, etc. For 
the hrst term in Eq. 04.15p . terms in H that have any positive powers of r will not contribute 
by the analysis of Sec. 13.1.11 For the second term we need to keep terms in H up to order r^, 
because the derivatives will reduce the order by 2. 

The symmetrical combination H{t, + , t), will lead to something whose Fourier trans¬ 

form does not decline rapidly for positive so that if this alone were put into Eq. fl4.4p the 
integral over ^ would not converge. But each term in + H{t',t) will combine with a 

term coming from iE{x,x') to give something whose Fourier transform does decline rapidly. 
We will work order by order in r. 

3.2.1 General computation of E 

We will need the Green’s functions for the background potential, including only hrst order in 
V, so we write 

G = + ■ ■ ■ . (3.25) 


4 Vi 


a! 


ff. 


(3.24) 
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The equation of motion is 


(□ + f/ (x))G{x, x') = {x — x'). 

Using x') = 5^^\x,x') and keeping only first-order terms we have 

{x, x') = —V(x)G^^^ (x, x '), 


(3.26) 


(3.27) 


so 


G^^\x,x') = - / d^x"G^^\x,x")V{x")G^^\x'',x'). 


For t > t" > t' we have for the retarded Green’s function, 




” |x"-x'| 2 ) = 


1 5{t'' - t' - 
An 


(3.28) 


(3.29) 


So we can write 


GSJ^(a;,a;') =-^ / - il„ '7,,, —^U(t",x"). (3.30) 


5(t"-F- |x"-x'|), 


Stt^ 


X" — X' 


Integrating over the second delta function we find t" = t' + |x" — x'|. Again considering purely 
temporal separation and defining Y' = x" — x' and Y' = we find 


- / dG / dC'C^ 


GYitY) = 


d(r2-2rC")..,., 


Stt^ 


C" 


U(F + C",x' + C"G), 


(3.31) 


where J dfl denotes integration over solid angle, and G varies over all unit vectors. We can 
integrate over to get (" = r/2 and 




(3.32) 


If we dehne a 4-vector G = (0, G) we can write 

41) +/ 




(3.33) 


The advanced Green’s functions are the same with t and t' reversed. Since E is the advanced 
minus the retarded function, we have 


= / d^V{x+'^-^n)sgnT. 


(3.34) 


3.2.2 Terms of order r ^ 

We now compute the various Hj, Hj, and Ej, starting with terms that go as a~^ or r“^. These 
terms are exactly what one would have for flat space without potential. Equation fl2.4p gives 


H_i{x,x') 


1 

An^a+{x, x') 


1 


47r2(r^ — Y) ’ 


(3.35) 


17 













where 

C = x-x' 

and 

C = IC|. 

Similarly, the advanced minns retarded Green’s fnnction to this order is 


E_i{x,x') = Ga{x,x') - Gr{x,x') 


- C) - + C) 


so 

1 1 , _5(r + C)-5(r-C) 

ri-e c 

where 

F(r+) = limFfr + ie). 

Taking the e —)■ 0 limit in r+ and r_ gives the formnla 


t') = lim —- 
^ c^o87r2 




. K't + C) - 

ITT -:- 


C) 


so 

as discnssed in Ref [15]. 


(3.36) 

(3.37) 

(3.38) 

(3.39) 

(3.40) 

(3.41) 

(3.42) 


3.2.3 Terms with no powers of r 

We can hnd the Hadamard coefficients from the recnrsion relations, Eq. fl2.13p . fl2.14p . To hnd 
the zeroth order of the Hadamard series we need only vq. For flat space, a^a = —‘^Vabi.x" — x')^ 
and □cr = —8. Pntting these in Eq. (12.1311 we have 

(x" - a;')%,a + Wo = ^ , (3.43) 


Now let x" = x' + X{x — x') to integrate along the geodesic going from x' to x. We observe 
that 


So Eq. (13.431) gives 


dvo{x", x') 
d\ 

dvo{x", x' 


A 


= {X- x'Yvo^a{.x'\x') . 


/ // /X 

+ no(a; ,x) = 


or 


d\ ' ^ 4 ’ 

d{\vo{x\x')) _ V{x") 
dX “ 4 ’ 


from which we immediately find 




(3.44) 

(3.45) 

(3.46) 

(3.47) 
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Now we consider purely temporal separation so the background potential is 
Ar, x). We expand E in a Taylor series in r around 0 with t hxed. 

evaluated at {f + 

V{t' + Xt) = V{i) + r(A - -)Vt{t) + —(A - + ■ ■ ■ 

(3.48) 

We are calcnlatine the zeroth order so we keep onlv the hrst term of Eq. 113.481). and Eq. fl3.47p 

gives 


vo{t,t') = -V{t) + 0{T^) 

(3.49) 

and thus 


Ho{x,x') - r!). 

(3.50) 

and 


Ho{x,x') + Ho{x',x) - ^^^V{x)ln\T\. 

(3.51) 

We can expand V around x, 


v{x + ^n) = v{x) + v^^\x + ^n ), 

(3.52) 

where is the remainder of the Taylor series 


r r 

V^^\x + -n) = V{x + -n)-V{x)= / drVi{x + rn)QX 
2 2 Jo ' 

(3.53) 

Then from Eq. (13.341). 


Eo{x,x') = ■^E(F)sgnr 

OTT 

(3.54a) 

Ro{x,x') - j dnV^^\x+^^^n)sgnT. 

(3.54b) 

Using 

2 In r + TTZsgnr = In (—r^), 

(3.55) 

we combine Eqs. (13.511). (l3.54aD to End 


Ho{t,t') - r!). 

(3.56) 

Combining all terms through order 0 gives 


H(o){t,t') = + HQ{t,t') + ^iRo{t,t'). 

(3.57) 


3.2.4 Terms of order 

Now we compute the terms of order in H and E. First we need Uq at this order, so we use 
Eq. fl3.48p in Eq. fl3.47p . The V^t term in Eq. fl3.48p does not contribute, because it is odd in 
A — 1/2, and the others give 

vq{x,x') = ^V{x) +^-• (3-58) 

4 9o 
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Next we need to know vi, but since vi is multiplied by in H, we need only the r- 
independent term Vi{x,x). From Eq. fl2.14p . 

(□ + 1/(a;))no(a;, x') + 2?7“V,a(x, x')afi{x, x') + vi{x, x')'^xcr{x, x') = 0 . (3.59) 

We neglect the V{x)vq term because it is second order in V. At x = x', = 0, so 

ni(x,x) = ^ lim □a,no(x,x'). (3.60) 

O x'^x 


Using Eq. fl4.62p we hnd 

^xVo{x,x') = - [ dXU^V{x'+ X{x — x')) = - [ dA A^(nU)(x'+ A(x — x')), 
4 Jo 4 Jo 

and Eq. fl3.60l) gives 


ni(x,x) = —□U(x) 
96 


We also need to know wi, but again only at coincidence. Reference 

wi{x, x) = 

Combining the second term of Eq. fl3.58p with Eqs. fl3.62p . fl3.63|) gives 


gives 


3 1 

■-vAx.x) =- nvix). 

2 ^ ’ 64 ^ ^ 


HAtA') = 


1287r2 


^Uii(x)ln(-r!) + ^□U(x) 


(3.61) 

(3.62) 

(3.63) 


(3.64) 


Then Hi{x',x) is given by symmetry, so 
idi(x, x') + HAx', x) = 


r 


647r2 


-Uii(x) In |r| + -□U(x) 


(3.65) 


The calculation of Ei is similar to that of Eq, but now we have to include more terms in 
the Taylor expansion of V around x. So we expand 


U(x + Jfl) = U(x) + Wax^t + + V^^\x + , 

Z Z o Z 

where the remainder of the Taylor series is 


(3.66) 


1/(3) (x + -12) = - f ^ dr VijAx + rl2) ^ n^Xim^dr . (3.67) 

«J 0 

Since f dflQ^ = 0 and f dflfl^fl^ = (47r/3)W, Eq. (I3.34p gives 

Ei(x,x') = Y^l/ii(x)r^sgnr, (3.68a) 

Ri(x,x') = j dQV^^\x. (3.68b) 

Again using Eq. fl3.55p . we combine Eqs. fl3.65p . fl3.68ap to get 

= ^\n{-Tl)Vii + ^nV{x) . (3.69) 

Combining all terms through order 1 gives 

%)(t,T) = H-AtA') + HoitA') + HAtA') + ImAA') ■ (3.70) 
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3.3 The 


Using Eqs. fl3.9p . fl4.15p . we need to compute 


where 


= g{t)g{t') 


TT 




Using Eqs. fl3.42p . fl3.54bp . fl3.56p . fl3.57p . fl3.68bl fl3.69p . fl3.70p we can write this 


F{t, f) = g{t)g{t') ^ fi{t, t') , 


i=l 


with 


27r^rl 


u = 

128^2 [K..(u + v,m\ 


h = 

2567r2 

j dnvl 

vW{x+'-^n) 

i sgn r 

k = 

1 

647r2 

j dndl 

(x + ■^^)* Sgll T 


3.4 The Fourier transform 


(3.71) 

(3.72) 

(3.73) 

(3.74a) 

(3.74b) 

(3.74c) 

(3.74d) 

(3.74e) 

(3.74f) 


We want to calculate the quantum inequality bound B, given by Eq. fl3.8p . We can write it 


where 




i=l 


B, = 




Br = 
Bs = 
using Eqs. fl2.20p . 


poo Jp poo poo 

/ — dt dt'g{t)g{t')fi{t,t')e^‘ 

*0 ^ J—oo J—oo 

(jf ^ T ^ T 

/ ~ dtg{t--)g{t+-)fi{t,T)e-^^^,i = l 

Iq ^ J—oo J—oo ^ ^ 

/ oo 1 poo 

dtg\t)Q{t)= I 


..6 


7687r2 


dtg\t)DV{t) 


/ dtg^{t)CVii{t), 

' —OO 

and fl3.63p . 


(3.75) 


(3.76a) 

(3.76b) 

(3.76c) 
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3.4.1 The singular terms 

For i = 1,2, 3, /j consists of a singular function of r times a function of t (or a constant). So 
we will separate the singular part by writing 

= gi{t)si{T). (3.77) 

Then we dehne 

/ OO — ^ 

digi{^g{i--)g{t+-), (3.78) 

B,= [ drG',(r)s,(r)e-'«T (3.79) 

Jo J—OO 

This is a Fourier transform of a product, so we can write it as a convolution. The Gi are all 
real, even functions, and thus their Fourier transforms are also, and we have 

1 roo roo 

Now if we change the order of integrals we can perform another change of variables g = ^ 
so we have 

1 noo p—oo ^ POO 

= j dvGiiv)siiC) = ^ j dgGi{g)hi{g), (3.81) 

where ^ 

hi{g) = [ dCsiiC ). (3.82) 

The arguments of Ref. |T5] show that the integrals over ^ in Eq. fl3.79l) and g in Eq. fl3.82p 
converge. 

We now calculate the Fourier transforms in turn, starting with Bi. We have 


3 

».(«) = 2^^ 

(3.83a) 


(3.83b) 

The Fourier transform of Si is [20] 


^''i(C) = fc^0(C), 

(3.84) 

so 

hi{g) = • 

(3.85) 

From Eq. (I3.8ip we have 

1 r°° 

®‘=24j„ 

(3.86) 

Using 7'(0 = *^/(0> we get 

(3.87) 
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The function Gi is even, so its Fourier transform is also even and we can extend the integral 

1 r . —. . 1 


Bi = 


dSvr 


For Gi we have 


Gi(r) = 


27r2 


^hGr(r7) = -Gr(0). 


dtg{t-^)g{t+^) , 


and taking the derivatives and integrating by parts gives 

/*oo 


Bi = 




dtg"{t}^ , 


reproducing a result of Ref. 
For B 2 we have 




S2[T) = 

This calculation is the same as before except the Fourier transform of S 2 is 

S2(C) = -27rC0(C). 

So we have 


where 


G2{t) = 


87r2 


After taking the derivatives 


Bo = 


327r2 


B 2 = ^G'm, 


dtV{t)g{t-'^)g{t+^] 


dtV{t}[g{t}g"{t} -g'it}'^] 


For i ?3 we have 

S3(r) = ln(-r!). 

In the appendix, we find the Fourier transform of S 3 as a distribution, 

noo 

S 3 [/] = dvr / dk f'{k) In \k\ - 47r7/(0). 

Jo 

From Eq. f|3.82|) . we can write 

/ CX> 

dC soiO^iv - C), 

-00 

which is given by Eq. fl3.97l) with f{() = — C)y so 

poo 

hs{ri) = —At: / d(^ 5{r] — () In \(\ — ATT'yQ{ri) = —ATTQ{ri){\nT] + 7) . 


(3.88) 


(3.89) 


(3.90) 


(3.91a) 

(3.91b) 

(3.92) 

(3.93) 

(3.94) 

(3.95) 

(3.96) 

(3.97) 

(3.98) 

(3.99) 
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Then Eq. fl3.8ip gives 


B. = - 


TT 


c/? 7G3(?7) (Inr^ + q) =-/ d7]G‘i{r]){\n\7]\ + '-f), 


TT 


(3.100) 


since G 3 is even. The integral is just the distribution w of Eq. fIB.lOp applied to (^ 3 , which is 
by definition ta[G 3 ], so Eq. fIB.lip gives 


B. = 


dr In |r| sgn r , 


(3.101) 


with 


G 3 (r) = 


OOtt^ 




so 


B. = 


OOvr^ 


drln|r|sgnr / dtV^uil^git-'^)g'{t+'^). 


(3.102) 

(3.103) 


3.4.2 The non-singular terms 

For i = 4, 5, 6, fi is not singular at r = 0. We include everything in 

/ OO 

di fiir,i)gii--)g(t+-), (3.104) 

B,= r drF,(r)e-'«" = H ^F,(-0 = - T d^0(OF,(-O . (3.105) 

Jo J—00 J0 J—00 

The integral is the distribution 0 applied to Fi{—^)), which is the Fourier transform of 0 
applied to Fi^—r). The Fourier transform of the 0 function acts on a function / as [20] 

©[/] = ^ ’ (3.106) 

where P denotes principal value, so 

B, = ~P J dr (^i"*(r)^ + F,( 0 ). (3.107) 

The first of the non-singular terms is a constant: does not depend on r. Thus F 4 is even 

in r, and only the second term of Eq. fl3.107p contributes, giving 

1 roo 

B, = F4(0) = j dtg(tf IKwffl + K«(i)] ■ (3.108) 

The functions /s and /g are odd in r and vanish as r ^ 0, so in these cases only the 
hrst term in Eq. fl3.107p contributes and the principal value symbol is not needed. Equa¬ 
tions fl3.74ell3.104113.107p give 

1 7°° 1 _ _ r _ r 7 Irl 

^5 = y dr-J dtg{t--)g{t+-) j dVLVlV^^\x+^-jVL)sgnT (3.109) 
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and Eqs. fl3.74fjl . fl3.104p and fl3.107p give 


B. = 


647r3 


dr^ f dtg{t-^)g{t+^) j dVtdl 


{x + sgn T 


(3.110) 


Here we can integrate by parts twice, giving 
1 


B. = 


647r^ 


dr 


' —OO J —OO 


didl -g{t-^)g{t+^) 

Til 


dll (x + ^11) sgnr . (3.111) 


From Eqs. (13.531) . fl3.67p . goes as r and as for small r, so the r integrals converge. 


3.5 The Quantum Inequality 


Now we can collect all the terms of B from Eqs. fl3.76bp . fl3.76cp . fl3.90p . fl3.95p . fl3.103p . fl3.108p . 
(13.1091) and (I3.11ip . Since H 7 is made of the same qnantities as H 4 , we merge these together. 
We hnd 

1 


B = 


where 


IOtt^ 


h + -I2 - -I3 - -il + - —iX 

2 ^ 6 ^ 8 ^ lOvr ® dvr ® 


-L 


V 

7 ) 


(3.112) 


/i = 


- 

Ir, - 


- 


-*4 


4^ 


- 


rV 


' —OO 

poo 


' —OO 

poo 


' —OO 

poo 


’ —OO 

poo 


dtg"{tf 

diV{l)[g{t)g''{t) - ^'(^^] 


7~ T 

diVii(t)g(t- -)g'(t + -) 

dtgit)'^ 

0 0 


dr- 


dtg{t-^)g{t+I dnVlV^^\x + 


’ —OO 

poo 


T 


J —OO 
'OO 


dr 


dtdi 




dfl {x H-fl) sgn 


. T 


= C dtg{t}‘^Vii{i). 


(3.113a) 

(3.113b) 

(3.113c) 

(3.113d) 

(3.113e) 

(3.113f) 

(3.113g) 


In the case with no potential, only Ji remains, reproducing a result of Fewster and Eveson [9]. 

In Eq. (I3.113cl) . In |r| really means ln(|r|//), where I is the arbitrary length discussed in 
Ch. [21 The choice of a different length changes Eqs. (I3.113cp and (]3.113dp in compensating 
ways so that B is unchanged. 

Equations (II.2113. 112113.113p give a quantum inequality useful when the potential V is known 
and so the integrals in Eqs. (13.1131) can be done. If we only know that V and its derivatives 
are restricted by the bounds of Eq. (I3.10p . then we can restrict the magnitude of each term of 
Eq. (I3.1T2P and add those magnitudes. We start with 

/ OO ^00 

/ dt[g{t}\g"{t}\+g'it}'^] . (3.114) 

■OO J —OO 
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The cases of , Ij, and are similar. For and Ig , it is nsefnl to take explicit forms for 
the Taylor series remainders. From Eq. fl3.53p . we see that 




< 


I dSJivn^iiin-i < J = 9’^I^I'CL- 

(3.115) 


Similarly from Eq. fl3.67p we have 


dnv^^\x+'-^n) 


< 


< 


|r| 


T 


48 


dn\Vijk\\n^nm^\ 

■ vZ^Y f 

ijk 


(3.116) 


27r+ 1 


rlV'" , 

' I * max ’ 


We can then perform the derivatives in Eq. fl3.113fl) and take the absolnte valne of each resnlting 
term separately. 


We define 


J2 

Js 

Ja 

Jr> 

Je 

Jr 


dt [^(f)|/(f)| +g'{tY] 

/ oo 

dt'\g'{t')\g{t)\\n\t' - t\\ 

-OO 

poo 

dtg{tf 

/ OO 

dt'g{t)g{t') 

■OO 

poo 

dt I dt'\g'{t')\g{t)\t' -t\ 
dt I dt' [g{t)\g"{t') \ + g'{t)g'{t')] {t' - tf 


’ —OO 

poo 


' —oo 
poo 


' —oo t/ —oo 


' —oo J —oo 


' —oo J —oo 


and find 


'2 I — 


'3 1 ^ 

'4 I — 
'5 1 — 


fdnax'-^2 

< ^^maxJs 

< 97rE"Lj; 


tVi 


‘6 I — 


< 


max"^ 5 

— ^ (4-^5 + 4J6 + Jr) 


Wl 

-'7 I 


< 3|C|E"^J4 


(3.117a) 

(3.117b) 

(3.117c) 

(3.117d) 

(3.117e) 

(3.117f) 

(3.118a) 

(3.118b) 

(3.118c) 

(3.118d) 

(3.118e) 

(3.118f) 


Thns we have 

[ dtg{t)^{Tr)4t,0) > 


IGvr^ 


dl + -ldnax'42 + K 


// 

max 


+E' 

I y r 


/// 

max 


ll7r + 1 

IGvr 


J5 + 


2 

277 + 1 
6477 


11 


o'^s T ( 777 T 4877 \C\ ) J 4 


24 

(4J6 + Jr) 


.(3.119) 
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3.5.1 An example for a specific sampling function 

An example of the quantum inequality with a specific sampling function g is the following. 
Consider a Gaussian sampling function 

a(t) = . (3,120) 

where to is a positive number with the dimensions of t. Then the integrals of Eqs. fl3.117p . 
calculated numerically, become 


Ji = 3.75to^ 

J 2 

= 3.15to^ 


J 3 = 2.70to 

J 4 

= 1.25to 

(3.121) 

J 5 = 3.14f2 

Je 

= 3.57f^ 


J 7 = 3.58t ^, 





so the right hand side of Eq. fl3.119p becomes 

- {3-75 + + (1.63 + 591.25|C|)V;;„tS + . (3.122) 

3.6 Discussion of the result 

In this chapter we have demonstrated a quantum inequality for a flat spacetime with a back¬ 
ground potential, considered as a first-order correction, using a general inequality derived by 
Fewster and Smith, which we presented in Ch. [2] . We calculated the necessary terms from 
the Hadamard series and the antisymmetric part of the two-point function to get H. Next we 
Fourier transformed the terms, which are, as expected, free of divergences, to derive a bound 
for a given background potential. We then calculated the maximum values of these terms to 
give a bound that applies to any potential whose value and first three derivatives are bounded. 

To show the meaning of this result, in the last section we presented an example for a specific 
sampling function. By studying the result we can see the meaning of the right hand side of our 
quantum inequality. The first term of the bound goes as where to is the sampling time, 
and agrees with the quantum inequality with no potential [9]. The rest of the terms show the 
effects of the potential to first order. These corrections will be small, provided that 

VtraJl < 1 (3.123a) 

« 1 (3.123b) 

K'Vi; « 1. (3.123c) 

Equation fl3.123ap says that the potential is small when its effect over the distance to is 
considered. Given Eq. fl3.123ap . Eqs. fl3.123bD and fl3.123cp say, essentially, that the distance 
over which V varies is large compared to to, so that each additional derivative introduces a 
factor less than Unlike the flat spacetime case the bound does not go to zero when the 
sampling time to oo so we cannot obtain the Averaged Weak Energy Gondition (AWEG). 

Finally, it is interesting to note the relation of this result to the case of a spacetime with 
bounded curvature. Since the Hadamard coefficients in that case are components of the Rie- 
mann tensor and its derivatives, we expect that the bound will be the flat space term plus 
correction terms that depend on the maximum values of the curvature and its derivatives, just 
as in our case they depend on the the potential and its derivatives. We will demonstrate that 
this hypothesis is true in the next chapter. 
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Chapter 4 


Quantum Inequality in spacetimes 
with small curvature 


In this chapter we present a derivation of a timelike-projected quantum inequality in spacetimes 
with curveture. First we consider a massless, minimally-coupled scalar field with the usual 
classical stress-energy tensor, 


= (4.1) 

Let 7 be any timelike geodesic parametrized by proper time t, and let g{t) be any any smooth, 
positive, compactly-supported sampling function. 

Let’s construct Fermi normal coordinates j27] in the usual way: We let the vector eo{t) be 
the unit tangent to the geodesic 7 , and construct a tetrad by choosing arbitrary normalized 
vectors 6^(0),i = 1,2,3, orthogonal to eo(0) and to each other, and define {ei(t)} by parallel 
transport along 7 . The point with coordinates (a:°, x^, x^) is found by traveling unit distance 

along the geodesic given by x^ei{0) from the point 7(0). 

We work only in first order in the curvature and its derivatives, but don’t otherwise assume 
that it is small. We assume that the components of the Ricci tensor in any Fermi coordinate 
system, and their derivatives, are bounded. 


I I ^ Rmax 


\R. 


ab,cd\ — 


< R'l 


\R. 


I < R'" 

ab,cde\ -'T'max 


These lead to bounds on the Ricci scalar and its derivatives. 


|R| <4R„ 


R,cd\ < 


\R,cde\ < 4Rmax > 


(4.2a) 


(4.2b) 


since we are working in four dimensions. 

Eqs. fl4.2|) are intended as universal bounds which hold without regard to the specihc 
choice of Fermi coordinate system above. We will not need a bound on the hrst derivative. 
The reason that we bound the Ricci tensor and not the Riemann tensor is that, as we will 
prove, the additional terms of the quantum inequality do not depend on any other components 
of the Riemann tensor. 

Thus we can write Eq. fl2.23l) in our case as 




(4.3) 






where 


-B — / + f dtg^it) (q — 2aR^ii — —{Rtt,tt + Rii,tt — ^Rtt,n + Rii,jj) ) ) (4-4) 


and F denotes the Fonrier transform in both argnments according to Eq. 02.241) . 


(4.5) 


4.1 Simplification of 

The operator of Eq. 02.3p for a massless held can be written 

1 


ry-iSplit _ 

Ft' - 2 


dtdf + ^ didi! 


2=1 


(4.6) 


To simplify it, we will dehne the following operator, 

3 




( 4 . 7 ) 


2 = 1 


which in hat space would be the derivative with respect to the center point. Then Eqs. 04. 6 p 
and 04.7p give 


psplit 

■tt' 


8A' + 5 (V? - - v=,) 


— ^ [v? + Da; — + Da,/ — d^i + 2dtdti^ , 


(4.8) 


where Da; and Da;' denote the D’Alembertian operator with respect to x and x'. Because we are 
using Fermi coordinates and are on the generating geodesic, the D’Alembertian and Laplacian 
operators have the same form with respect to Fermi coordinates as they do in hat space. Then 
using 

(4.9) 


= 7 + dl] , 


we can write 




□ a;Tf + ^x'H + Viiiif 


dlH. 


( 4 . 10 ) 


Consider the hrst term. The function H{x, x') obeys the equation of motior^l in x and so does 
E{x,x'). Thus 

n,H =^0,H{x',x). (4.11) 

As we discussed in Sec. 13.1.21 we have 


H{x', x) = H{x, x'Y + {wj{x\ x) — Wj{x, x'))aYx, x '), 


( 4 . 12 ) 


^In general the sums in Eq. (12.41) do not converge and we should work only to some finite order in a. In 
that case H{x,x') obeys the equation of motion to that order. 
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Since is real, O^H{x,x')* = 0, and we have 

= ^n^'^{wj{x',x)-Wj{x,x'))a^{x,x'). (4.13) 

j 

In the coincidence limit Eq. fl4.13p vanishes. There is no j = 0 term becanse tco = 0. In the 
j = 1 term, we have a, which vanishes at coincidence nnless both derivatives of the □ are 
applied to it, in which case wi cancel each other, and for j > 1, even vanishes. 

The second term in brackets in Eq. (I4.10p gives 

= ^n^,'^{wj{x,x') -Wj{x',x))a^{x,x'). (4.14) 

j 


Adding together Eqs. fl4.13p and f|4.14l) . we get something which is smooth, symmetric in x 
and x', and vanishes in the coincidence limit. Following the analysis of Sec. 13.1.21 snch a term 
makes no contribntion to Eq. (14.4p and for our purposes we can take 


y split ^ 



(4.15) 


4.2 General computation of E 

The function E is the advanced minus the retarded Green’s function, 

E{x^ x') = Ga{x^ x') — Gji{x, x '), (4.16) 


and iE is the imaginary, antisymmetric part of the two-point function. The Green’s functions 
satisfy 

aG(xy) = ^EE^, (4,17) 

Following Poisson, et ah [30] and adjusting for different sign and normalization conventions. 


G{x, x') = — {2U{x, x')6{a) -I- V (x, a;')0(—a)) , 
4:71 

where U{x,x') = A^^‘^{x,x') and V{x,x') are smooth biscalars. 

For points y null separated from x', V is called V [30] and satisfies 


(4.18) 


+ 


-□a+ 2 


V = -□[/, 


(4.19) 


with all derivatives with respect to y. Now V is first order in the curvature, so we will do 
the rest of the calculation as though we were in flat space. Under this approximation, we will 
neglect coefficients which depend on the curvature, and also evaluate curvature components at 
locations that would be relevant if we were in flat space. The distance between these locations 
and the proper locations is hrst order in the curvature, so the overall inaccuracy will always 
be second order in the curvature and its derivatives. 
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Thus we use (T’“ = 


2{y — x'Y and Dcr = —8 in Eq. fl4.19p to get 


[y-xrYa{y) + V{y) = ]^m{y). 


(4.20) 


Now suppose we want to compute V at some point x". We need to integrate along the geodesic 
going from x' to x”. So let y = x' + \{x" — x') and observe that 

= x'^YY + ^^y'l = Kx"-xrv, + v(y) = {y-x'rv^ + v{y) = \aU{y), ( 4 . 21 ) 

SO ^ 

V{x\x') = ]^j^ dXUU{y). (4.22) 

The function V obeys [30] 

n^V(x,x') = 0. (4.23) 

Let us consider points x and x' on the geodesic 7 , which in the flat-space approximation means 
they are separated only in time, let x = (x-l-x')/2, and establish a spherical coordinate system 
{r,6,(j)) with origin at the common spatial position of x and x'. Now V{x,x') can be found 
in terms of V and its derivatives evaluated at the time t (the time component of x) using 
Kirchhoff’s formula, 


V {x, x') = — / dQ 


V(x",x') + ^§;y(x",x') + 


(4.24) 


where the derivatives act on the first argument of V, f dQ means to integrate over all spatial 
unit vectors 12 , and we now set 

x" = x + (r/2)12 (4.25) 

with the 4-vector 12 given by 12 with zero time component. 

Now define null coordinates u = t + r and v = t — r. Then x" has u = t, v = 0. The 
derivative d/du can be written {d/dt + d/dr)/2 and so 

V{x,x') = ^ j [ME((V2)12i,a:')]„=^ , (4-26) 

where 12i is 12 with unit time component. From Eq. fl4.22p . 

7/ 1 /■“ 

uV{-n,x') = -J^ du\aU){{u'/2)n,,x'), (4.27) 

with the D’Alembertian applied to the first argument, and so 

V(x,x') = — [ dil0^ifU{x”,x'). (4.28) 

8vr J 

We are only interested in the first order of curvature, so we can expand U, which is just 
the square root of the Van Vleck determinant, to first order. From Ref. |33j . 

A^/^(x,x') = 1 — - / ds{l — s)sRab{sx + {1 — s)x'){x — x'Y{x — x'Y + 0{RY , (4.29) 

2 Jo 


31 










so in the case at hand we can use 

U{x\ x') = x') = l-- / ds{l - s)sRah{y)X‘^X^ (4.30) 

2 Jo 

where y = sx” = (sm", sn", 6*", 0") is a point between 0 and x'\ and the tangent vector X = 
dy/ds. We are interested in 0). To bring the □ inside the integral, we dehne Y = 

sX = (sm", sn", 0, 0), and 

nu{x",0) = -l [ ds{l-s)sn,4Rab{y)X^X^] = ~ [ ds{l - s)sny[R,,{y)Y^Y^]. (4.31) 
^ Jo ^ Jo 

For the rest of this section, all occurrences of m, n, 6*, 0, and derivatives with respect to these 
variables will refer to these components of y or Y. 

Now we expand the D’Alembertian in Eq. fl4.28p . in terms of an angular part, 

_2 4 9 f . 

^ {v-uysinem 

and a radial and temporal part, which we can write in terms of derivatives in u and v, 




+ 


02 


{v — uY sin 


(4.32) 


02 


0 0 


dvdu u — V \du dv 
The angular part vanishes on integration, leaving 


V{x^x') = —— f dQ f dss{l — s) 


dud,, - 


u — V 


{du - 0 ,;) 


{Rab{y)Y^Y^) 


Outside the derivatives, we can take n = 0 and change variables to m = sr, giving 

1 


V{x,x') = 


dvrr^ 

1 


dQ J du{T - u) [ududu - 0« + 0^] {Rab{y)Y‘^Yy 
dVt [ du{T - u)du[{udy - l){Rab{y)Y^Yy]. 


dvrr^ _ 

We can integrate by parts with no surface contribution, giving 

1 


V{x,x') = 


1 

47J-T-3 


dnj du{l - udu){Rab{y)Y‘^Yy 

dfl / duu [ uRuu^ui^y') 2,Ruu(^y') T Ruu(^y')^ 


Now 


Rab — Gab “ (l/2)5'abO 5 

where Gab is the Einstein tensor and G its trace. Thus 
1 


V{x,x') = 




dVt / duu [—uGuu,v{y) — “^Guviy) Y y/‘^)G{y) + Gy 


Now dehne a vector held Qa{y) = Gab{y)Y^. Then 

Qa;c = GabAy)Y' + Gab{y)Y^-,c. 


(4.33) 

(4.34) 

(4.35) 

(4.36) 

(4.37) 

(4.38) 

(4.39) 

(4.40) 
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We write the covariant derivative only becanse we are working in nnll-spherical coordinates, 
rather than becanse of spacetime cnrvatnre, which we are ignoring becanse we already have 
hrst order qnantities. 

Since the covariant divergence of G vanishes, 

g'^'^Qa;c = g'^'^Gab{y)Y’^;c. (4.41) 

In Cartesian coordinates, = y^, and = 5^, which means that (in any coordinate system). 

g^^Qa;c = G. (4.42) 


Explicit expansion gives 


g Qa\c — ‘^{QvjuYQujv)' 


U — V 


iQ u Qv 


{v — uy 


Id 1 

(sin 6 ^( 50 ) + 


sin 6 do 


sin6'^ 


, (4.43) 


bnt the angnlar terms vanish on integration. Now we expand the derivatives in u and v and 
set n = 0 , giving 


Qv,U ^Guv,U Y Guv 

Qu,V ^Guu,V 4 “ Guv 1 


so 


J dQ (2uGuv,u + ‘2uGuu,v + 8Guv — 4:Guu) — J dVL G. 

Snbstitnting Eq. fl4.45p into Eq. fl4.39|) . we hnd 

V{x,x') = J J [i^GuvAy) + ‘^Guviy) - Guuiy)] 

and integration by parts yields 


v{x,x') = — dn 


Guv{x") — ^ duu^ (Guviy) + Guuiy)) 

Jo 


Now 


dVL j duu^ [Guviy) + Guuiy)) = n duu"^ [Guiy) + Gtriy)) 


(4.44a) 

(4.44b) 

(4.45) 

(4.46) 


(4.47) 


dn duu^G^\y)-G^^iy)) (4.48) 


which is 4 times the total flnx of crossing inward throngh the light cone. Since this qnantity 
is conserved, = 0 , we can integrate instead over a ball at constant time i, giving 


dfl 


r/2 


dr r^G“(x + rfl) = — j dVL j ds 5 ^( 7 “ (x + s(r/2)f2) 


4 


(4.49) 

















where x" = a; + s(r/2)r2, and 


Gr{x,x') = A^^‘^{x,x')^j^ + 


271 327r2 


dn{ - [Gu{x'') - Grrix")] 


ds s'^Gttix'!.) 


(4.51) 


r^/ /N A 1/2/ - 1^ - + |X - X'l) 

E{x,x) = A'^'^{x,x)— -^ 


4 ?! X — x' 


(4.52) 


327r2 


1 


- [Gttix") - Grr{,x'')] - ds s^Gtt{,x”) [> sgn r . 


4.3 Computation of H 

We now need to compnte H{x, x') and apply First we consider the term in H{x, x') that 

has no dependence on the curvature. It has the same form as it would in flat space as shown 
in Ref [T5] and Ch. [3] 

H-i{x,x') = H_i{x,x') = ^ --. (4.53) 

A7r^a+[x, x'j 

In Sec. 14.41 we will apply the fully general from Eq. fl4.15p with defined in Eq. fl4.7p 
to H_i{x, x'). 

All the remaining terms that we need are hrst order in the curvature, so for these it is 
sufficient to take as the flat-space Laplacian with respect to the center point, x. For this 
we only need to compute H at positions given by time coordinates t and t' but the same spatial 
position. 

As we discussed in Sec 13.1.11 we only need to keep terms in H with powers of r up to r^, 
but we need E exactly. The terms from H alone give a function whose Fourier transform does 
not decline fast enough for positive ^ for the integral in Eq. fl4.4p to converge. Thus we extract 
the leading order terms from iE and combine these with the terms from H. This combination 
gives a result that has the appropriate behavior after the Fourier transform. 


4.3.1 Terms with no powers of r 

First we want to calculate the zeroth order of the Hadamard series. The Hadamard coefficients 
are given by Eqs. fl2.13112.14p for a massless field. To find the zeroth order of the Hadamard 
series we need only vo{x,x'), which we find by integrating Eq. fl2.13p along the geodesic from 
x' to X. Since we are computing a hrst-order quantity, we can work in flat space by letting 
y' = x' + \{x — x') and using the hrst-order formulas na = —8 and cr’“ = —2{y' — x'Y- From 
Eq. fl2T3P . we have 

(//' - x'fvo^a + Vo = ^□A^/2(//', x '), (4.54) 

and thus 

1 G 

vq{x,x') = - / (iA(nA^/^)(a;'-I-A(a; — x'), a:'). (4.55) 

4 Jo 

by the same analysis as Eq. fl4.22p . 
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Using the expansion for from Eq. fl4.29p gives 


"1 




vo{x,x') = -- d\ ds{l - s)say,[R^^{sy'+ {1 - s)x'){y'- x'y{y'- x'f] 
° Jo Jo 

1 rl rl 


dX / ds{l — s)s 


(As)^(ni?ab)(x' + sX{x — x')){x — x')°‘{x — x') 


/\b 


+2XsRb{x' + sA(a; — x)){x — x')^ + 2R{x' + sA(a; — x')) 


We can combine the s and A integrals by dehning a new variable a = sX 




a a 
~X^^ 


d(x - af 


(4.56) 


(4.57) 


(4.58) 


Then, changing a to s, we hnd 


Vo{x,x') = / cis(l-s) 

16 Jo 


^ {URg^f^){x' + s{x — x')){x — x'y{x — x')’’ 

+2sR^b{x' + s{x — x')){x — x')’’ + 2R{x' + s{x — x')) 


or when the two points are on the geodesic, 


vo{t,t') =-—j ds{l-sf 


(Di^tt) {x' + sr)r^ + 4:sr]‘^Rct,d{x' + st)t + 2R{x' + st) 


In the second term we nse the contracted Bianchi identity, ri^^Rct,d = R^t/‘Ji giving 


2 [ ds{l — sYsrRtix'+ st) = 2 [ ds{l — s)'^s—R{x'+ st) 
Jo ’ Jo ds 

= —2 / (is(l — s)(l — 3s)i?(a:'+ sr), 

Jo 


so the hnal expression for no is 


vo{t,t') =-—j ds{l-s) 


s^(l - s)aRtt{x + (s - l/2)r)r^ + 4si?(x + (s - l/2)r) 


(4.59) 


(4.60) 


(4.61) 


(4.62) 


To calculate Hq we only need the zeroth order in r from no, so the hrst term does not contribute. 
In the second term, we make a Taylor series expansion, 

R{x + (s - l/2)r) = R{x) + i?,t(h)r(s - 1/2) + ^R^ttix)T‘^{s - 1/2)^ + O(r^), (4.63) 

but only the hrst term is relevant here. Thus 


vo{tR') = -T / ds{l - s)sR{x) = -—R{x). 


24 


(4.64) 
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We also need to expand the Van Vleck determinant appearing in the Hadamard series. From 
Eq. draD, 

t') = l- ^Ru{x)t^ - + OA ). (4.65) 

Keeping the first order term from Eq. fl4.65p and using Eq. fl4.64p . we have 


Ho{x,x') = 


1 


487r2 


Rtt{x) - -R{x) In(-r^) 


Now we can add the Ho{x',x) which is the same except that t and t' interchange 

Ho{x, x') + Hq{x\ x) = iRtt(x) - R(x) In |r|] . 


(4.66) 


(4.67) 


Next we must include E from Eq. fl4.52p . We can expand the components of the Einstein 
tensor around x, 

GaA”) = Gab[x) + {x") , (4.68) 

where is the remainder of the Taylor series 


g1;) (x") = GAx”) - Gab{x) = I dr G,A^ + • 


r/2 


(4.69) 


To hnd Eq, we put the first term of Eq. fl4.68|) into the second term of Eq. fl4.52p . We use 
Grr = Gij^VVt^ and j = (47r/3)W and hnd 


Eq{x,x') = —l-Gtt{,x)--Gii{x) 


1 f 1 

)7r 

^G(x)sgnr = -^i?(x)sgnr. 
487r 487r 


ds s^Gtt{x) \ sgn T 


(4.70) 


For the remainder term Rq, we put the second term of Eq. fl4.68p into the second term of 
Eq. g32]), 


Rq{x,x') = 


327r2 


dVti- 


1 r 


G^u A') - ( x")l - ds s^gH'^ (x'J) Ugn r . (4.71) 


Using 

2 In |r| + vrisgnr = In (—A), 
we combine Eqs. fl4.67|) and fl4.70|) to hnd 


(4.72) 


HQ{t,t') = 


1 


487r2 


Rtt{x) - ^R{x) ln(-r!) 


(4.73) 


Combining all terms through order 0 gives 


1 


= i7_i(f,F) + + -iRo{t,t'). 


(4.74) 
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4.3.2 Terms of order 

Now we compute the terms of order in H and E. To find vq at this order we take Eqs. (I4.62p 
and fl4.63p and include terms through second order in r. The hrst-order term vanishes, leaving 


Vq{x,x') = ds{l - s)[s'^{l - s)nRtt{x) 




+2s(s - l/2f Ru{x)y + ... (4.75) 

1 


480 


(^Rtt{x) + -R^ttix)^ + ... . 


Next we need vi but since it is multiplied by in H we need only the r independent term. 
From Eq. fl2.14p 

□uo + 2ni^aCr’“ + niDcr = 0 , (4.76) 

At X = x', (T’“ = 0 so 


ni(x, x) = - lim n^vo{x,x'). 


x—>x' 


(4.77) 


Using Eq. fl4.59p in Eq. fl4.77p . the only terms that survive in the coincidence limit are those 
that have no powers of x — x' after differentiation, so 


1 /■! 1 
ni(x, x) = —— / ds{l — sYs‘^OR(x) =———OR{x). 


16 


480 


(4.78) 


Equations fl4.62p . fl4.75p and fl4.78p agree with Ref. [6] if we note that their expansions are 
around x instead of x. 

The Wi at coincidence is given by Ref. 


(4.79) 


3 1 

u;i(x,x) = --ni(x,x) = • 


Combining Eqs. fl4.75l) . fl4.78p . and fl4.79p . and the fourth order term from the Van Vleck 
determinant of Eq. fl4.65l) . and keeping in mind that a = —when both points are on the 
geodesic, we hnd 


i7i(x,x') = 


6407r^ 


-Rtt,tt{^) - - o In (-r!) 


(4.80) 


Then Hi{x',x) is given by symmetry so 


i7i(x, x') + i7i(x', x) = 


1607r2 


6 


Rtt,tt{.x) - -□R(x) - - ( □Rn(x) + -R,uix) I In |r 


r 


(4.81) 

The calculation of Ei is similar to Eq, but now we have to include more terms to the Taylor 
expansion, 

43) ('4_g2) 


Gabix") = Gabix) + ^GabA^)n^ + (^) + G]^’A") ^ 

Z o 


where the remainder of the Taylor series is 


gSU") = 2 


/X \ ^ • j. 

drGab,ijk{x + rfl) y- - rj . 


(4.83) 
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Now we put Eq. fl4.82p into Eq. fl4.52p . and again use Grr = . The first term of 

Eq. fl4.82p gives Eq, which we computed before, and the second term gives nothing, be¬ 
cause J = f = 0. Using J and f = 

(47r/15)(5*-^^^^ -|- -|- the third term gives 


Ei{x,x') = 


1927r 




ds s'^Gtt,ii{x) 


T sgnr 


(4.84) 


3207r 


^Giijj{x) + ^Gij^ij{x) - ^Gu,ii{x) sgnr . 


1 

-( 

2 


Using the conservation of the Einstein tensor, 0 = ri°‘°Gia,b = Gu^t — Gijj and 0 = y]°‘°Gta,b = 
Gtt,t — Git^i we can write 

So 


^ij,ij^x') Giiii{x) 


Now Gab = Rab - (l/2)i?, SO 


G,, = (3/2)i?ii - (l/2)i?,, 

Gtt = (l/2)i?ft -|- {l/2)Rii. 


Putting these in Eq. (I4.86p gives 


Ei{x,x') = (^Riijji^) + ^Rtt,ttix) + ^Rii^ttix)^ r^sgnr] 


QbOvr 

1 

9607r 


DRn{x) + ^R,tt{x) ) r^sgnr] . 


The fourth term of Eq. fl4.82p gives the remainder 


(4.85) 


Ei{x,x') = + Gtt,tt{x) - ^Gu^ii{x)^ r^sgnr. (4.86) 


(4.87a) 

(4.87b) 


(4.88) 


i?l(x, x') = 


327r2 


dfl 


2 L 


Gif (i") - G<f (i" 


dss^Gif«)lsgnr. (4.89) 


To calculate Hi, we combine Eqs. fl4.81l) and fl4.88p and use Eq. fl4.72p to get 


Hi{x,x') = 


6407r^ 


-Rtt,tt{x) - -□7?(x) - - ( URn{x) + -R,tt{x) In(-r^) 


All terms through order 1 are then given by 


(4.90) 


Hm{t,t') = H-i{t,t') + HQ{t,t') + Hi{t,t') + -tRi{t,t '). 


(4.91) 
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4.4 The 

We can easily take the derivatives of Hq and Hi nsing Eq. fld.lSp . because they are already 
hrst order in R. However in the case of the term we have to proceed more carefully. 

From Eqs. fl4.7p and fl4.53p we have 



1 A / a" 2 --^ V M 

47r2 dx^ dx'"^ d{x''^Y) V'^+y 


1 av av \ 

47r2(T+ dx'^dx'^ ’ 


(4.92) 


where we used da/dx'^ = da/dx'^ = 0 when the two points are on the geodesic. From [6j, after 
we shift the Taylor series so that the Riemann tensor is evaluated at x, we have 


d{x^Y 

d'^a 

d{x'^Y 

av 

dx^dx'^ 


2 1 1 

= -2riii - -Riut{x)T'^ - -RititA^)r^ - -Ritit,ttr^ + 0 {t^) 

2 1 1 

= -2r]ii - -Ritit{x)T^ + -Ritit,t{^)T^ - -Ritit,ttT^ + 0{t^) 

2r]ii - ^Riutix)T^ - ^Ritit,ttT'^ + O(t^) . 


From Eqs. fl4.92p and fl4.93p . and using Rmt = —Rtt we have 


= 




2 3 

+ —Rtt,tt{x) 


(4.93a) 

(4.93b) 

(4.93c) 


(4.94) 


From Eqs. fl4.5p and fld.lSp . we need to compute 



TT 




(4.95) 


where 


= 9{t)g{t') 




(4.96) 


In the hrst term in brackets it is sufficient to use i7(o)(t,T); because higher order terms in H 
are smooth, even in r, and vanish at coincidence, and so they do not contribute, as discussed 
in Sec. 13.1.11 In the second term, two powers of r are removed by differentiation, so we need 

Using Eqs. (I4.53p . (14.711) . (I4.73p . (14.741) . (I4.89p . (14.901) . (14.911) and (14.941) we can combine 
all terms in F to write 

6 


F{t, t') = g{t)g{t') ^ fi{t, t') , (4.97) 

i=l 


39 

















































with 


f2 


/s 

h 


fb 


k 


3 

2n‘^T‘t 


1 

487rV2 


[Rii{x) - 7Rtt{x)] 


(4.98a) 

(4.98b) 


1 


3847r2 

1 

3207r2 

1 


'll 3 

~^Rtt,tt{x') ~\~ ~Rii,tt(^x') Rtt,iiik) ~^Rii,jj^ 

43 __ 7 _ 1 


X) 


ln(-r!) 


2567r2 / ^ 2 L 


Ga^(k') - (x")l - / (x") I i sgn r 


-^ / d^dl , 

647r2 / ^ 2 L 


G^^{x")-GfXH- / dss^Gfix'kUsgnr 


(4.98c) 

(4.98d) 


(4.98e) 


(4.98f) 


4.5 The quantum inequality 

We want to calculate the quantum inequality bound B, given by Eq. fl4.4p . We can write it 

8 

B = J2Bi, (4.99) 

i=l 


where 


B^ 


By 

Bs 


OO /•CO 


— / dt I dt’g{t)g{t')fi{tR')R^G t) 


' —OO J —OO 
/•CO /•OO 


— / dr I dtg{t-'^)g{t+'^)fi{t,T)e = 1 __Q 


TT 


' —OO «/ —OO 


^OO -1 

2 / - ^ -L 


dtg {t)Q{t) = 


38407r2 


dtg^{t)DR{i) 


r-oo T 

2 -' „ /-X 0 


dtfii (t) 


2ci/?jj(x) ^(yRiiii{x) -\- Riiiiik) 2>Riiii{x') 


RRii,ji{k)) 


(4.100a) 

(4.100b) 


(4.100c) 


using Eqs. fl2.20p . fl2.22p . fl4.4p and fl4.79p . The hrst 6 terms have the same r dependence as the 
corresponding terms in Ch. [31 So the Fourier transform proceeds in the same way, except that 
instead of dependence on the potential and its derivatives, we have dependence on the Ricci 
tensor and its derivatives. After the Fourier transform, we see that R 4 and B-j have exactly 
the same form so we merge them in one term. Thus 


B = 


lOvr^ 


T \ ^ tR ^ I ^ rR I ^ ^ tR 

' + 12 2 24 + 240 ^ + lOvr ' 47 r ^ 
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-'7 ' 


(4.101) 
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where 
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(4.102b) 
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(4.102f) 


dt^ (t) 


2cid? jj(x) ^iyRtt,ttiyX') 4“ Rii^ttix') 3Rti^ii{x^ 




(4.102g) 


If we only know that the Ricci tensor and its derivatives are bounded, as in Eqs. fl4.2p . we can 
restrict the magnitude of each term of Eq. fl4.10ip . We start with the second term 


tR\ 


‘2 I — 


< 


dt \7Ru{x) - Rii{x) \ \g{t)g''{t) - g'{t}g'{i)\ 


< lOR, 


dt[g{t}\g"{r> \ + g'{t) ]. 


(4.103) 


Terms and are similar. For and Iq, we need bounds on the components of 

G. From Eq. fl4.87bp . \Gtt\ < 2,Rj^^^. Since Eq. fl4.2l) holds regardless of rotation, we can 
bound Grr at any given point by taking the x-axis to point in the radial direction. Then 
Grr = Gxx = (1/2)[Rxx “ Ryy “ Rzz + Rtt] and \Grr\ < 2,R^gj^. Taking derivatives of G just 
differentiates the corresponding components of R, which are also bounded. In particular, since 
there are 3 terms in V|, we have |V|Git|, |V?Grr| < QR' ^^^ . Using these results and Eq. fl4.69p 
for the remainder we have 


j davl 


- gHHx'') 



< Uydfijl [|v2G„,i(i)| + |V"G„,i(i)|] +ybss=|WG„,i(i)| 




(4.104) 


41 
































For Iff we use Eq. fl4.83|) for the remainder 


dn{- 


1 r 


+ / dss^Gfix'i) 


< 


\T\ 


48 


^ [\Grr,ijkix)\ + \Gtt,jk{x)\] + I ds s^\Gu,^Ad^)\ )■ Iff 


< M 
— -^max 


id,k 


3 D'" 


rhi?: 


After we bound all the terms and calculate the derivatives in we can dehne 
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7(27r + l) 
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<ax(4J5 + 4J6 +Jr) 


rR\ 


< (24|a| + ll|6|)7?LxJ4. 


Thus the hnal form of the inequality is 

f dTg(tf{TrUt,0)>-j^J 
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(4.105) 

(4.106a) 

(4.106b) 

(4.106c) 

(4.106d) 

(4.106e) 

(4.106f) 

(4.107a) 

(4.107b) 

(4.107c) 

(4.107d) 

(4.107e) 

(4.107f) 

(4.108) 


J 3 + ( ^ + IStt (24|a| + 11|7|) ) J 4 


1637r + 14^ 7(27r + l),^^ 

J 5 + ff 4J6 + J7) 


OOvr 


1927r 


Once we have a specihc sampling function g, we can compute the integrals of Eqs. fl4.106p 
to get a specihc bound. In the case of a Gaussian sampling function, 

g{t) = , (4.109) 
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we computed these integrals numerically in Sec. 13.5.11 Using those results, the right hand side 
of Eq. fld.lUSp becomes 


- {3.76 + + |1.71 + 197.9(24|o| + ll|6|)]fl"„io + 6.99fl'lx‘o} ■ (^-HO) 

The leading term is just the flat spacetime bound of Ref. |9] for g given by Eq. fl4.109p . The 
possibility of curvature weakens the bound by introducing additional terms, which have the 
same dependence on as in Ch. [31 with the Ricci tensor bounds in place of the bounds on 
the potential. 

4.6 Discussion of the result 

In this chapter, using the general quantum inequality of Fewster and Smith we presented in 
Ch. 121 we derived an inequality for a minimally-coupled quantum scalar field on spacetimes 
with small curvature. We calculated the necessary Hadamard series terms and the Green’s 
function for this problem to the first order in curvature. Combining these terms gives H and 
taking the Fourier transform gives a bound in terms of the Ricci tensor and its derivatives. 

If we know the spacetime explicitly, Eqs. fl4.3l) . fl4.10ip . and fl4.102|) give an explicit bound 
on the weighted average of the energy density along the geodesic. This bound depends on 
integrals of the Ricci tensor and its derivatives combined with the weighting function g. 

If we do not know the spacetime explicitly but know that the Ricci tensor and its hrst 
3 derivatives are bounded, Eqs. (14.1061) and fl4.1U8|) give a quantum inequality depending on 
the bounds and the weighting function. If we take a Gaussian weighting function, Eq. fld.llOp 
gives a bound depending on the Ricci tensor bounds and the width of the Gaussian, to. 

As expected, the result shows that the corrections due to curvature are small if the quanti¬ 
ties Rmax^O) -^max^O) -^max^o much less than 1. That will be true if the curvature is 

small when we consider its effect over a distance equal to the characteristic sampling time to 
(or equivalently if to is much smaller than any curvature radius), and if the scale of variation 
of the curvature is also small compared to to. 

In all bounds, there is unfortunately an ambiguity resulting from the unknown coefficients 
of local curvature terms in the gravitational Lagrangian. This ambiguity is parametrized by 
the quantities a and b. 

Ford and Roman [13] have argued that flat-space quantum inequalities can be applied in 
curved spacetime, so long as the radius of curvature is small as compared to the sampling 
time. The present chapter explicitly confirms this claim and calculates the magnitude of the 
deviation. The curvature must be small not only on the path where the quantum inequality is 
to be applied but also at any point that is in both the causal future of some point of this path 
and the causal past of another. All such points are included in the integrals in Eq. fl4.1U2ep 
and fl4.102fl) . 

Is is interesting to consider vacuum spacetimes, i.e., those whose Ricci tensor vanishes. 
These include, for example, the Schwarzschild and Kerr spacetimes, and those consisting only 
of gravitational waves. In such spacetimes, the flat-space quantum inequality will hold to 
first order without modihcation. There are, of course, second-order corrections. For the 
Schwarzschild spacetime, for example, these were calculated explicitly by Visser [331133]. 137]. 
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Chapter 5 

Average Null Energy Condition in a 
classical curved background 


In this chapter we present the proof of the achronal ANEC in spacetimes with curvature 
using a null-projected quantum inequality. It is structured as follows. First we state our 
assumptions and present the ANEC theorem we will prove. We begin the proof by constructing 
a parallelogram which can be understood as a congruence of null geodesic segments or of 
timelike paths. Then we apply the general inequality presented in Ch. |2] to the specihc case 
needed here, using results from Ch. 01 Finally we present the proof of the ANEC theorem 
using the quantum inequality. 

5.1 Assumptions 

5.1.1 Congruence of geodesics 

As in Ref. m, we will not be able to rule out ANEC violation on a single geodesic. However, 
a single geodesic would not lead to an exotic spacetime. It would be necessary to have ANEC 
violation along a finite congruence of geodesics in order to have a physical effect. 

So let us suppose that our spacetime contains a null geodesic 7 with tangent vector i and 
that there is a “tubular neighborhood” M' of 7 composed of a congruence of achronal null 
geodesics, defined as follows. Let p be a point of 7 , and let Mp be a normal neighborhood of 
p. Let n be a null vector at p, linearly independent of £, and let x and y be spacelike vectors 
perpendicular to v and i. Let q be any point in Mp such that p can be connected to g by a 
geodesic whose tangent vector is in the span of {v,x,y}. Let 7 (g) be the geodesic through g 
whose tangent vector is the vector £ parallel transported from p to g. If a neighborhood M' of 
7 is composed of all geodesics 7 (g) for some choice of p, Mp, v, x and y, we will say that M' 
is a tubular neighborhood of 7 . 

5.1.2 Coordinate system 

Given the above construction, we can define Fermi-like coordinates described in Appendix lAl 
on M' as follows. Without loss of generality we can take the vector v to be normalized so 
that Va£^ = 1, and x and y to be unit vectors. Then we have a pseudo-orthonormal tetrad 
at p given by E(u) = £, E(^) = v, E(^x) = x, and E{jj) = y. The point g = {u, v, x, y) in these 
coordinates is found as follows. Let g^^^ be found by traveling unit affine parameter from p 
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along the geodesic generated by +yE{jjy Then q is fonnd by traveling nnit affine 

parameter from gd) along the geodesic generated by uE^^). During this process the tetrad is 
parallel transported. All vectors and tensors will be described using this transported tetrad 
unless otherwise specihed. 

The points with u varying but other coordinates hxed form one of the null geodesics of the 
previous section. 


5.1.3 Curvature 

We suppose that the curvature inside M' obeys the null convergence condition, 

RabV’^V^ > 0 (5.1) 

for any null vector l/“. Equation fl5.ip holds whenever the curvature is generated by a “classical 
background” whose stress tensor obeys the NEC of Eq. fll.3l) . We will refer to this as a “classical 
background”, but the only way it need be classical is Eq. fll.3p . 

We would not expect any energy conditions to hold when the curvature is arbitrarily large, 
because then we would be in the regime of quantum gravity, so we will require that the 
curvature be bounded. In the coordinate system we described we require 


\R 


abed 


< Rn 


(5.2) 


and 


\R, 


abed,a \ 


< RL 


\R 


■abed,a 


/?! <r: 


// 

max’ 


I-Ro.6cd,a/37 1 R 


in 

max 


(5.3) 


in M', where the greek indices 13^'y ■ ■ ■ = v,x,y and i?max, -Rmax) -^max) R'max finite numbers 
but not necessarily small. Also we assume that the curvature is smooth. 


5.1.4 Causal structure 

We will also require that conditions outside M' do not affect the causal structure of the 
spacetime in M' 1110 

J+ (p, M) n M' = J+ (p, M') (5.4) 

for all p G M'. Otherwise the curvature outside M' may be arbitrary. 


5.1.5 Quantum field theory 

We consider a quantum scalar field in M. We will work entirely inside M', and there we require 
that the field be massless, free and minimally coupled. Outside M', however, we can allow 
different curvature coupling, interactions with other fields, and even boundary surfaces with 
specified boundary conditions. 

Because M may not be globally hyperbolic, it is not completely straightforward to specify 
what we mean by a quantum field theory on M. We will use the same strategy as Ref. m- 
Our results will hold for any quantum field theory on M that reduces to the usual quantum 
field theory on each globally hyperbolic subspacetime of M. The states of interest will be those 
that reduce to Hadamard states on each globally hyperbolic subspacetime, and we will refer 
to any such state as “Hadamard”. See Sec. II B of Ref. HH for further details. 

^This condition is equivalent to J“(p, M) n M' = J~{p, M') for all p G M'. 
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(|)= (u,v,0,0) 


(0,v,0,0) 



(0,0,0,0) 


Figure 5.1: Construction of the family of null geodesics <F using Fermi normal coordinates 

5.2 The theorem 

5.2.1 Stating the theorem 

Theorem 1. Let {M,g) be a (time-oriented) spacetime and let 7 be an achronal null geodesic 
on {M,g), and suppose that 7 is surrounded by a tubular neighborhood M' in the sense of 
Sec. 15.1.11 obeying the null convergence condition, Eq. fll.Sp . and that we have constructed 
coordinates by the procedure of Sec. 15.1.21 Suppose that the curvature in this coordinate 
system is smooth and obeys the bounds of Sec. I5.1.3( that the curvature in the system is 
localized, i.e., in the distant past and future the spacetime is flat, and that the causal structure 
of M' is not affected by conditions elsewhere in M, Eq. fl5.4p . 

Let uj he a state of the free minimally coupled quantum scalar held on M' obeying the 
conditions of Sec. 15.1.51 and let T be the renormalized expectation value of the stress-energy 
tensor in state u. 

Under these conditions, it is impossible for the ANEC integral, 

/ OO 

dXTa,ri\TiX)), (5.5) 

•CO 

to converge uniformly to negative values on all geodesics F(A) in M'. 

5.2.2 The parallelogram 

We will use the (u, v, x, y) coordinates of Sec. 15.1.21 Let r be a positive number small enough 
such that whenever |n|, |x|, \y\ < r, the point {0,v,x,y) is inside the normal neighborhood Np 
dehned in Sec. 15.1.11 Then the point (u, v, x, y) G M' for any u. 

Now consider the points 

$(«, n) = (u, n, 0, 0). (5.6) 

With V hxed and u varying, these are null geodesics in M'. (See Fig. 15.11 ) Write the ANEC 
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integral 

/ OO 

duTuu{^{u,v)). (5.7) 

-OO 

Suppose that, contrary to Theorem 1, Eq. fl5.7p converges uniformly to negative values for all 
|n| < r. We will prove that this leads to a contradiction. 

Since the convergence is uniform, A{v) is continuous. Then since A{v) < 0 for all |n| < r, 
we can choose a positive number vq < r and a negative number —A larger than all A{v) with 
V G (— no,no)- Then it is possible to hnd some number ui large enough that 



duTuu{,^{u,v)) < -y4/2 


(5.8) 


for any v G (—Vo, "^0) as long as 


M+(n) > Ml (5.9a) 

U-{v) < —Ml. (5.9b) 

As in Ref. m, we will dehne a series of parallelograms in the (m, m) plane, and derive 
a contradiction by integrating over each parallelogram in null and timelike directions. Each 
parallelogram will have the form 


V G (-Mo, Mo) (5.10a) 

M G (m_(m), m+(m)) , (5.10b) 


where m_(m), m+(m) are linear functions of v obeying Eqs. (15.Oh . On each parallelogram we will 
construct a weighted integral of Eq. (15. 8 h as follows. Let /(a) be a smooth function supported 
only within the interval (—1,1) and normalized 



(5.11) 


Then we can write 


rVQ 


dvf{v/vof 


i-u+iv) 


' —VO 


f U— ( r ) 


duTuu{^{u,v)) < -V0AI2. 


(5.12) 


We can construct this same parallelogram as follows. First choose a velocity V. Eventually 
we will take the limit R —?■ 1. Dehne the Doppler shift parameter 


Let a be some hxed number with 0 < a < 1/3 and then let 

to = . (5.14) 

As R —)■ 1, (5 —)■ cx) and to —)■ 0. 

Now dehne the set of points 

= . (5,15) 
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Figure 5.2: The parallelogram <F(m,u), v G (—Uo,Uo); 
t G (-to,^o), V e (-ho,ho) 


G (m_(u),M+( u)), or equivalently 


We will be interested in the paths given by *Fy (h, t) with r] hxed and t ranging from —to to to- 
In flat space, such paths would be timelike geodesic segments, parameterized by t and moving 
at velocity V with respect to the original coordinate frame. In our curved spacetime, this is 
nearly the case, as we will show below. Dehne 


ho = ui + to6/V2 

(5.16a) 

Vo = to/{V25) 

(5.16b) 

u±{v) = ±ho + S‘^v 

(5.16c) 


so that u± satishes Eqs. fl5.9p . Then the range of points given by Eq. fl5.6p with coordinate 
ranges specihed by Eqs. fl5.10p is the same as that given by Eq. fl5.15p with coordinate ranges 


— to < 
-ho < 


The parallelogram is shown in Fig. 15.21 
The Jacobian 


d{u, v) 


d{r],T) 


so Eq. fl5.12p becomes 


t < to 
h < ho 


1 



dTTuu{^v{v,t))f{'t/toY < 


—Ato/2. 


We will show that this is impossible. 


(5.17a) 

(5.17b) 


(5.18) 


(5.19) 


5.2.3 Transformation of the Riemann tensor 

Since we are taking 6 —)■ oo, components of R with more u’s than n’s diverge after the trans¬ 
formation. Components of R with fewer u’s than n’s go to zero and components with equal 
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numbers of u’s and n’s remain the same. We want the curvature to be bounded by i?max in the 
primed coordinate system, which will be true if all components of the Riemann tensor with 
more u’s than n’s are zero. We will now show that this is the case in our system. 

All points of interests are on achronal null geodesics, which thus must be free of conjugate 
points. Using Eq. fl5.ll) and proposition 4.4.5 of Ref. [22], each geodesic must violate the 
“generic condition”. That is to say, we must have 

i^i^^[aRb]cd[eif] = 0 (5.20) 

everywhere in M'. 

The only nonvanishing components of the metric in the tetrad basis are guv = Qvu = — 1 
and Qxx = Qyy = 1- The tangent vector £ has only one nonvanishing component £“ = 1, while 
the covector has only one nonvanishing component ^v = “I- Thus Eq. fl5.20p becomes 


^[aRb]uu[e^f]=^- (5.21) 

Let j, k, I, m and n denote indices chosen only fromU {x,y}. Choosing a = m, e = n, and 
b = f = V we hnd 

Rmuun = 0 (5.22) 

for all m and n. Thus 

Ruu = 0. (5.23) 


Equation fl5.23p also follows immediately from the fact that since Ruu cannot be negative, any 
positive Ruu would lead to conjugate points. 

If we apply the null convergence condition, Eq. (15.ip . to U = + eE(m) + (e^/2)E(„), 

where e <C 1, we get 

Ruu + ‘^Rmu^ + C(e^) > 0 . (5-24) 


Since Ruu 
have 

Since Rmu 


0 from Eq. fl5.23l) . in order to have Eq. fl5.24p hold for both signs of e, we must 


Rm.u 0 


-Rv 


T g^ Rjmku: 


(5.25) 


R 


umvu 


g^’'Rjmku ■ 


(5.26) 


Now we use the Bianchi identity. 


R 


luum;n 


R, 


lunu\m 


R, 


iumn\u 


= 0 . 


(5.27) 


From Eq. fl5.22l) . Riuum,n = 0. The correction to make the derivatives covariant involves terms 


of the forms RauunNnEi 


(a) 


and RiaurrNnE^\ Because of Eq. (15.221) . the only contribution 


to the hrst of these comes from a = v, which we can transform using Eq. fl5.26p . For the 
second, we observe that 0 = Vn{E^'""> ■ = 2VnEu\ so a = v does not 

contribute. Furthermore Riumn-,u = Riumn,u, because the u direction is the single hnal direction 
in the coordinate construction of Sec. 15.1.21 and so in this direction the tetrad vectors are just 
parallel transported. Thus we find 

dRlumn 


du 


= gE[R^mku^nEr + Rjlku^nEt^ - R^uku^mEr “ Rjlku^mE^^^] 
EiRlkum + Rlukm)^nE^^ + {Rlknu + Rlunk)'^niE^^ . 


(5.28) 


^This notation applies only in this subsection and not the rest of the thesis 
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Eq. fl5.28p is a first-order differential equation in the pair of independent Riemann tensor 
components Rxuxy and Ryuxy By assumption, the curvature and its derivative vanish in the 
distant past, and therefore the correct solution to these equations is 

Rlumn = 0 . (5.29) 

Eqs. fl5.26p and fl5.29p then give 

Rumvu — 0 • (5.30) 

Combining Eqs. fl5.22p . fl5.29p . and fl5.30p and their transformations under the usual Rie¬ 
mann tensor symmetries, we conclude that all components of the Riemann tensor with more 
u’s than n’s vanish as desired. 

5.2.4 Timelike paths 

The general quantum inequality of Ch. [2] we will use for this proof is applied on timelike paths. 
So we are going to show that the paths ^{r]+5t/ \/2, t/\/2(5) are indeed timelike. Differentiating 
Eq. fl5.15p . we find the components of the tangent vector k = d^y/dt = (l/\/2, l/\/2) in the 
Fermi coordinate basis. The squared length of p in terms of these components is gabk°‘k^. 
We showed in Appendix that gab = Pab + hab, where hab at some point X is a sum of a 
small number of terms (6 in the present case of 2-step Fermi coordinates) each of which is a 
coefficient no greater than 1 times an average of 

RabcdX^X^^ (5.31) 

over one of the geodesics used in the construction of the Fermi coordinate system. The sum¬ 
mations over d and c in Eq. (15.31 p are only over restricted sets of indices depending on the 
specific term under consideration. From Eqs. (I5.15p and (I5.16al) the points under consideration 
satisfy 

|m| < ui/S + -\/2tq (5.32a) 

\v\ < ro/\/2 (5.32b) 

X = y = 0. (5.32c) 

From Eq. fl5.14p . the hrst term in Eq. fl5.32al) decreases faster than the second, so we hnd 
that all components of X are O(to)- Using the fact that the components of the Riemann tensor 
are bounded we hnd 

hab = O(Rmax^o) (5.33) 

so 

gabk°'k^ = —1 + O(Rmaxto) • (5.34) 

Thus for sufficiently large 5, and thus small A; is timelike. 

5.2.5 Causal diamond 

The quantum inequality of Ch. [2]Kontou:2014tha is applied to timelike paths inside a globally 
hyperbolic region of the spacetime. So this region N must include the timelike path from 
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p = ^viVy “^o) to g = $^(77,^0), and to be globally hyperbolic it must include all points in 
both the future of p in the past of q, so we can let N be the “double cone” or “causal diamond”, 

N = J+{p)nJ-{q) . (5.35) 

We have shown that the curvature is small everywhere in the tube M', so we must show that 
N C M'. 

From the previous section, we have that the metric can be written as 

9ab Vab j (5.36) 

where hab consists of terms of the form Rabcd.X'^X^. The double cone in flat space obeys 

|a;|, I2/I, |^^| < To , (5.37) 

so the same is true at zeroth order in the Riemann tensor R. Thus at zeroth order, 

hab = 0{R^aA ), (5-38) 

and so at hrst order in R, 

|x|, \y\, |n| < ro(l + O(Rmaxto)) • (5-39) 

Since R <^r for large 5, we have 

|a:|, \y\,\v\ < r . (5.40) 

and N C M' as desired. 


5.3 The null-projected quantum inequality 


We can write the general quantum inequality of Eq. fl2.23p for w(t) = ^yijqR) for a specihc 
value of r] and the stress energy tensor contracted with null vector held £“ = n as 



dTg(tY(TZ)(^(t)) > -B, 


( 5 . 41 ) 


where 


B= 

Jo ^ 


dtg^{t) {2aR^uu + 


where we used that Quu = 0 so the term Q doesn’t contribute at all. Run - 
Sec. 15.2.31 and 

F denotes the Fourier transform in both arguments according to Eq. fl2.24l) . 


(5.42) 
0 according to 

(6.43) 


5.4 Calculation of 

uu' I-*-! 


To simplify the calculation we will evaluate the TFIiy^ in the coordinate system (f, a:, t/, z) where 
the timelike path w(t) points only in the t direction, z direction is perpendicular to it and x 
and y are the previously dehned ones. More specihcally t and z are 


5 + 5v 5 — 5v 

75 “’ 75 “ 


(5.44) 
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The new null coordinates u and v are defined by 

t — z 

and are connected with u and n, 


t + z 


M = 5 , V = 6v. 

The operator can be written as 


(5.45) 


(5.46) 


(5.47) 


If we define ( = z — z' and u as the u coordinate of x, the center point between x and x' we 
have 



(5.48) 


5.4.1 Derivatives on H_i 

For the derivatives of i7_i it is simpler to use Eq. fl5.47p . We have 

du'duH-i = ^ (~) • 

In flat spacetime it is straightforward to apply the derivatives to i7_i. However in curved 
spacetime, there will be corrections hrst order in the Riemann tensor to both a and its deriva¬ 
tives. 

We are considering a path w whose tangent vector is constant in the coordinate system 
described in Sec. 15.1.21 The length of this path can be written 

I dx^ dx^ I - 

s{x,x^) = j 9ab{w{X)) — —= J dXy/gab{x'')Ax^Ax’^. (5.50) 

where Ax = x — x' and x" = x' + A Ax since dx°‘/dX is a constant. 

Now a is the negative squared length of the geodesic connecting x' to x. This geodesic 
might be slightly different from the path w. However, the deviation results from the Christoffel 
symbol T^^, which is 0{R). Thus the distance between the two paths is also 0{R), and the 
difference in the metric between the two paths is thus 0{R?‘). Similarly, the difference in length 
in the same metric due to the different path between the same two points is 0{R^). All these 
effects can be neglected, and so we take a = — 

Now using Eq. flA.26p of Appendix |A] we can write the hrst-order correction to the metric, 

9ab Vab T ^ab T ^ba j (5.51) 

where Fab is given by Eq. flA.28p of the same Appendix because the hrst step for x = j/ = 0 is 
in the v direction and the second in the u direction. By the symmetries of the Riemann tensor 
the only non-zero component is 

Faii{x")= [ dK{l-K)Riiavu{fix"^,x''^)x''^x''^ , (5.52) 

Jo 
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where we took into account the different sign conventions. Putting this in Eq. flS.SOp gives 


s(x,x') = / dX\/2Ax^Ax^ + 2FjjyAx'^Ax^ 


d\\/2 ( VAx^Ax^ + -Fiji,{Ax'^)^^^{Ax'^) . 




So to first order in the curvature, 


a{x,x') =—s{x,x') =—T^ + C —2 


dXFjjfjAx'" Ax'" 


We define the zeroth order a, 

and the hrst order, 
r(l) 


a^»>(x,I') = -r2+C^ 


a^^\x,x') = —2 dXFyyAx'"Ax" 


'0 


r-l /•! 


= -2 / dX dK{l-K)Ryuyu{Kx"'^,x''")x''^x"'^Ax"Ax^ 

Jo Jo 

= -2 f dX [ dy{(i-y)Ryuyu{y,x’'")Ax"Ax'", 

Jo Jo 

where we dehned ^ = x"'" and changed variables io y = k£. Now to hrst order. 


a 


( 1 ) 


cr+ a 


(0) 


(a 


( 0))2 ’ 


and the derivatives, 
d d 


dx'^ dx''^ J \a^ J 


C=o 


4 12 

=-1- 

^4 ^ ^6 ^5 


2\/2 / (1) (i)\ 1 




U ,U' 


o-;-,/-ro¬ 


ll) 


r. 


4 ,uu' 


Now we can take the derivatives of a , 


aJ'’ = -2 1 dXX^ I dy{i - y)Ryyyu{y,x"")Ax"Ax 


d 


A A 


di 


= -2 1 dXX I dyRuyuv{y,x'''")Ax'"Ax'" 


Similarly, 


= -2 dX{l-X)-^ dy{i-y)Ryyyu{y,x"^)Ax"Ax^ 
= -2 [ dX{l-X) [ dyRuyuy{y,x"")Ax"Ax". 


For the two derivatives of 


= -2 / i^A(l - X)XRuyuy{x''^ x"^)Ax"Ax" . 


(5.53) 


(5.54) 

(5.55) 

(5.56) 


(5.57) 


(5.58) 


(5.59) 


(5.60) 


(5.61) 
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(5.62) 


Now we can assume purely temporal separation, so Ax^ = Ax'" = t/\/2 and 


x" =—it"+ z,t"-z), 

^/2 

where z = [z + z')/2 and t" = t' + Xr. Then the derivatives of TT-i are 




C-2 


»1 


4-12/ dX dy{£ - y)Ryuyu{y,x”") 
Jo Jo 

nl 


-2\f2 I dX{l-2X) I dyRu^ui,{y,x"")T 




>0 


>0 


+ / dX{l-X)XRi,,^,{x"^x""y 


Let us dehne the locations = {k,x",x'") and 


< = + z),t" - z). 


Then Eq. (15.63^ can be written 




^-2 




dX 


12 / dK{l-K){x""yR,uvu{x'') 


+2\/2(l - 2A) / dKx”"Ruijuv{x'y)T - (1 - X)XRuyuv{x'') 


T 


The derivatives of iL_i can thus be written 

1/1 


'inr J- 


-2 


+ yi{t, r) ] + ^y 2 {t, r) + ^ysit, r) 


Tl \TT^ J r 

where the |//s are smooth functions of the curvature, 

2/i(/'?')= [ dXYyt”) y2{i,T)= f dX{l - 2X)Y2{t") yyt,T) = f dX{l 
Jo Jo Jo 

with 


yyn = 

Y2{t") = 


27r2 


27r2 


f dK{l - K){t" + zYRyuijuix'y), 
0 

d/t(t T z^Ryuvui^zY )) 


1 


ryt") = -_i?-,-,(a;") , 


where x” and x" are dehned in terms of t” by Eqs. fl5.62p and fl5.64p . 


(5.63) 


(5.64) 

(5.65) 

)■ 

(5.66) 

A ) Ay 3 ( t "), 

(5.67) 

(5.68a) 

(5.68b) 

(5.68c) 
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5.4.2 Derivatives with respect to r and u 


In Ch. m we calculated but for points separated only in time. Let us use coordinates 
(T, Z, X, Y) to denote a coordinate system where the coordinates of x and x' differ only in T. 
Ref. [26] gives 


where 


Hio){T,r) 

Ho{T,r) 


^_i(T, T') + Ho{T, T') + Hi{T, T') + T '), 

T') + Ho{T, T') + 1*Ro(T, T') , 


47r2(T-T'-ie)2 ’ 


487r2 

(T - Tf 

6407r2 


RTT{x)--R{x)\n{-{T-T'-tef) 


-Rtt,tt{x) — 


- URii{x) + -RM^) In (-(T - r - leY) 


(5.69) 

(5.70) 

(5.71a) 

(5.71b) 

(5.71c) 


The order-0 remainder term is 

1 


Ro(T,T') = 


327r2 


(El 


1 
2 L 


d^(x") - G‘Z(x") 

1 


'<(1) ( v/f\ 


(5.72) 


- / sgn(r-r'). 


where J (El means to integrate over solid angle with unit 3-vectors hi, the 4-vector 12 = (0,12), 
the subscript R means the radial direction, and we define X” = x + (1/2)|T — T'|12, X” = 
X + (s/2)|T — T'|12, and 


G%{X") = Gab{X") - Gab{x) = 


,\T-T'\/2 


dr Gab,i{x -1- rl2)12'^ 


(5.73) 


The order-1 remainder term is 

Ri(T,r) = 

^ ^ 327r2 


dl2 


1 
2 L 


G®,(A") - G™ (A'") 

1 


7(3) / TA//\ 


- / *s"GW(A;) Ugn(r-r'), 


(5.74) 

( 3 ) 

where G\^ is the remainder after subtracting the second-order Taylor series. We can write 


G® (A") = 


r\T-T'\/2 


drGAB,ijK{x + rl2) 


( 


T-T' 


V 2 


r . 


(5.75) 


When we apply the r and it derivatives from Eq. (I5.48p . we can take {T^Z^X^Y) = 
(t,z,x,y) and calculate dlHo, d^Ho, d^Hi, d^Ro, and d^Ri. Applying it derivatives to Hq 
gives 


d^Hn = 


487r2 


Rtt,Uu{^x') 2 ( ^—) 


(5.76) 
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For the derivatives with respect to r we have 



1 


R(x ), 


(5.77) 


and 


dlH, = 


3207r2 


+ ^Rit(x) ) (3 + In(-r):)) 


(6.78) 


in the r —)• 0 limit. 

Applying u derivatives to Rq gives 


d^.Rf) = 


1 r /•hl/2 

- I dVt I drdl 

327r2 J io 


— Grr,i{x'")] 


- J ds s^Gtt,i{x'”) sgn r , 


(5.79) 


where x'” = x + rfl and x'" = x + srfl. 

Now we have to take the second derivative of Ri with respect to r, which is T — T' in this 
case. This appears in three places: the argument of sgn in Eq. fl5.74p . the limit of integration 
in Eq. 05.751) . and the term in parentheses in Eq. 05.75p . When we differentiate the sgn, we get 
S{t) and S'{t). but since ~ r^, there are enough powers of r to cancel the S or d', so this 
gives no contribution. When we differentiate the limit of integration, the term in parentheses 
in Eq. 05.751) vanishes immediately. The one remaining possibility gives 



1 

1287r2 



[Gtt,ijk(x"') — Grr,ijk(x"')] 



dss Giiij}^{Xg 


sgnr. 


(5.80) 


5.4.3 Derivatives with respect to 

To differentiate with respect to C, we must consider the possibility that x and x' are not purely 
temporally separated. We will suppose that the separation is only in the t and z directions 
and construct new coordinates (T, Z) using a Lorentz transformation that leaves x unchanged 
and maps the interval (T — T', 0) in the new coordinates to (r, Q in the old coordinates. Then 


T -T' = sgnrv^r^ - (2 ^ 

and the transformation from (T, Z) to {t, z) is given by 


A = 


r C 

sgnr a/ r2 — \ C x 


with the X and y coordinates unchanged. Then 


r 


= A 


T-T' 

0 


(5.81) 


(5.82) 
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Now let M be some tensor appearing in The components in the new coordinate 

system are given in terms of those in the old by 

Mabc... = AaA^A^c • • • Mabc... ( 5 . 84 ) 

We would like to differentiate such an object with respect to ( and then set C = 0- The only 
place ( can appear is in the Lorentz transformation matrix, where we see 

O Aa I _ .^-1 / ^ ^ ^ /'K oka 


and similarly, 


^CA^ 


5?a:^ 


(5.85) 


(5.86) 


To simplify notation, we will dehne P and Q to be the matrices on the right hand sides. 
Reinstating x and y, 


P = 


Q = 


0 10 0 
10 0 0 
0 0 0 0 
0 0 0 0 

10 0 0 
0 10 0 
0 0 0 0 
0 0 0 0 


(5.87) 


(5.88) 


Now we can write the derivative of Mabc... as 


d(^MABC... — ■ ■ ■ )Afabc... 

C=o 


(5.89) 


= {d^AMh • • • + SlidcKWc ••• + ■■■ 1 

L J C=o 

= -{P16l,6^c • • ■ + ^PbS^C ••• + •• ■)Mabc... = -ptBC...Mabc... 

r_ N-/ r_ 

n 

where p'abc a rank-n matrix of O’s and I’s. With two derivatives, we have 


Q^Mabc... 


<=0 

= (92A“)5|5^--- + 5^(9^X)5^--- + ... 

+2(9cA^)(9cA^)5^ • • • + 2(9cA^)5|(9cA^c) ••• + ••• 1 Mabc... 


(5.90) 


T_ ^ 


+ P^Pl6^c--- + P16^sPc--- + ---)Mabc... 


a^b TDC 


(n—l)n 


= ^qTBC...Mabc... 
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where q'^BC ^ rank-n matrix of nonnegative integers. 

There are also places where T — T' appears explicitly in Hi. We can differentiate it using 
Eq. dmi), 


ac(T-T') 

dl{T-T') 


= 0 , 


C=o 


C=o 


= —r 


.-1 


(5.91a) 

(5.91b) 


Now we apply the operators 9? and d^dc to Hn, Hi, and H. First we apply one C derivative^ 

to Eq. (EUa usmg Eq. JS3a, 


I d(^Hii 


C=o> 


487r^r^ 


PttRabix), 


and two ( derivatives using Eqs. (I5.9np and (15. 91 ah . 

1 


d^^Ho 


C=o 




[quRabix) + R{x)] . 


(5.92) 


(5.93) 


Then we apply one ( derivative to Hi, 

1 


dr I 




19207r2 


Ptm Rab,cd{x) - Pu DR^h{x) + -Pu R,ab{x) (111 (-r_) + 2) 


and two ( derivatives to Hi, 
1 


(5.94) 


d^^Hi 


2=0 


6407r^ 


^QuuRab,cd{x) - ‘^RttA^) + OR{x) -^{q^aR^^{x) 
+^(lttR,ab{x) ) In (-r!) + ^ ( DRii{x) + ^R,tt{x) ) (1 + In (-r!)) 


(5.95) 

Finally we have to apply the derivatives to the remainder Ri. We can apply the ( derivatives 
in two places, the Lorentz transformations and Since the three terms are very similar 

we will apply the derivatives to one of them 


I dnePj^ix + Y) 


C=o 


= / dSJ I + Y)d^Y‘ 


(5.96) 


C=o, 


where we dehned = (1/2)|T — T'|A“r2^. Then using Eqs. fl5.89p and fl5.91ap . we hnd that 
that = (l/2)p“r2*sgnr and taking into account the properties of Taylor expansions, 




dY 




= Ggl(x"), 


(5.97) 


C=o 


3 The Lorentz transformation technique we use here is not quite sufficient to determine the singularity 
structure of the distribution at coincidence. Instead we can use Eq. (14.291) to compute the non-logarithmic 
term in Hq for arbitrary x and x', which is then —Rab{x){x — x')°‘{x — x')^/ (487r^cr+). Differentiating this term 
gives Eq. (15.921) and explains the presence of t_ instead of r in the denominator. The first term of Eq. (15.931) 
arises similarly. 
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(2) 

where is the remainder of the Taylor expansion of Gab,c after subtracting the first-order 
Taylor series. 

Thus Eq. fl5.96p becomes 




dQGP^ix + Y) 


C=o 


Using G^^^ from Eq. fl5.75p and 


d^l 



sgn T 


(5.98) 




drGab,ijc{x + rhi) 


Eq. fl5.98p becomes 




dVlG^^iX + Y)\^=o = 



, (5.99) 


+ ^PiGtt,ija{x'") sgnr 


-rj(5.100) 


We could simplify further by using the explicit values of the p matrices, but our strategy here 
is to show that all terms are bounded by some constants without computing the constants 
explicitly, since the actual constant values will not matter to the proof. 

Applying the r derivative gives 


dr 



dVlG^^ix + Y) 

\ r /•hl/2 

]= dn dr 


o 

II 


PttGab,ijkix'' 


+ ^PiGtt,ija{x" 




(5.101) 


We do not have to differentiate sgnr here, because the rest of the term is O(r^) and so a term 
involving 6{t) would not contribute. 

The same procedure can be applied to all three terms. Terms involving X'J will get an 
extra power of s each time G is differentiated. The final result is 


dr I dcRi{T,r) 
1 


(5.102) 


C=o, 


'hl/2 


327r2 


dfl 


dr 


^iPu-pTr)Gab,jk{x'")- / dsS^pt^Gab.jkK) 


Y{Gu,^Ja{x"') - Grr,ija{x'”)) - / 


Sgnr. 


For two ( derivatives we can apply both on the Lorentz transforms, both on the Einstein 
tensor or one on each, 

w d^ 


dl I dnGP^ix + Y) 


C=o 


j dn\ 

(<g_ 

^r2 

d 

^(3) 


^Ptt 

d 


dY^dY^ 


G^\x + Y)d^Y^d^Y^ 


C=o 

(5.103) 


r dY^ 


C=o 


C=o 
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Using Eqs. fl5.91bp and fl5.90p . 




C=o 


2r 2r 2r * ’ 


with hi = qj — 6l, while 




= 0 , 


C=o 


since qj = 0 and hi* = 0. Using properties of the Taylor series as before, we can write 

^2 




( 1 ) 


(5.104) 


(5.105) 


(5.106) 


so Eq. fl5.103p becomes 


[ dnG^Tix + Y) 

= 

[ dn 

J 

C=0 

J 


^ab 1 

T 




+j^(2pfeG® (x") + G'S,{^")h’ja 

Using as in Eq. fl5.73p and and from Eqs. fl5.99p and fl5.75l) this becomes 

>|t|/2 r /i \ 2 

= I dn I 

C=o 


(5.107) 


/ dnG^^^(x + U) 


/■|T|/2 r/i \ ^ 

dfl dr ( - - — j qt^GabM^"') 


+jPfp‘G„,.,.(x"') + (j - jj; 


+h\Gtt,ijk{x'‘ 




(5.108) 


For all three terms 

dlRiiT.T') 

1 


(5.109) 


C=o 


’M/2 


327r2 


dVt 


dr 


>0 


1 r 

2 ~ H 


2 (%t-C)GaM.;^(^'")- / dss\^,^Gat,jk{x'n 


'0 


~h~P'iP^j ni.Gtt,kab{x'") — Grr,kab{x'")) — f dss'^Gtt,kab{,x"J) 

Jo 


1 

'4^ 


r 


4 2|r| 


P-(Ptt -P“p)Ga6,,fcc(a:"') + ^-K{Gu,l,k{x''') - Grr,lA^"')) 


dss\2p<lptGab,kc{<') + KGu,m{<')) 


sgn r. 


5.5 The Fourier transform 

Eqs. (15T6D . ([526]), (jEnj)., (M), dEHOj), dESSj), dESSj), dSJl, (1^ . (I5T02D and 

(15.1091) include all the terms. To perform the Fourier transform we expand 
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according to Eqs. fl5.47p and fl5.48p and separate the terms by their r dependence, 




-2 


d^d^,H_^ + -{dlHo + -iR, 


~o(^T + + 29 t-5^) (Ho + Hi + -iRi 


= ^ ^ \+ yiitH)] + \y2{t,T) + +y3{t,T)) 


T \TT 




r 


+ In (-r!)c 2 (^ + coii) + C 4 (t, r) 


(5.110) 


where Ci, C 2 , and C 3 are smooth and have no r dependence and C 4 is odd, C\ and bounded. As 
mentioned in Sec. 15.4. 1[ the functions yi depend on r but are smooth. Explicit expressions for 
the Ci are 


Cl = 


C2 = 


C 3 = 


‘ ^ - fi(x) + H ^ ) R,,(x) 


1 


19207r2 

1 

960^ 


5/2, 


(5.111a) 

(x) + OR^,j(x) + - R^ab{x)^ (5.111b) 


5R 


■tt.uu I 


(S) - i (piu + 


2 / 


+p“ □i2a6(x) + i2,a6(a;) 


'hl/2 


(5.111c) 


C 4 = 


2567r2 


dr j dn{^l\^[Gtt,i{x''') - Grr,i{x''')]- ds s^Gu,i{x"') 


^ [Gu,ijkix"') - Grr,ijkix'")] dsS^Gu,ijkix7) 


4^2 

+ ( 1 - 7 ^ 


^(p?,^-pfr)Gab,M^"')- I dsS^pt,^GabMO 


2^ab/ 


+ ^-^{GuM^"') - Grr,jka{^'")) " / dsS^p^GuM^') 

^ .In 


,1 r 

+2 — 


T 


-{q?i - qfr)Gab,jk{x'") - / dss^qt^Gab,A<) 


1 


a^b 


-{GuMbW) - Grr,kab{x'")) - / dsS^GuMbiO 


T 


P,(P?? - pfr)GAikc{r") + pi‘AG,u,t(x"') - G„ijt{x"')) 


dss\2plptG^jtc(x":) + h‘fiut,k«)) ^a’a” )SJ-sgtir. 


J'l\ I 


k \ni , 


(5.111d) 


We now put the terms of Eq. fl5.110p into Eq. fl5.42p . and Fourier transform them, following 
the procedure of Sec. 13.41 , to obtain the bound B in the form 


b = s-^J2^^- 


(5.112) 


i=0 
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The first term in Eq. (15.1101) is l/(7r^rl), and we proceed exactly as Sec. 13.41 except for the 
different numerical coefficient, to obtain 


Bo = 


1 

247r2 



(5.113) 


Putting only Eq. fl5.113p into Eq. fl5.112p gives the result for flat space. Fewster and Eveson 
[9] found a result of the same form, but they considered Tu instead of Tuu, so the multiplying 
constant is different. Fewster and Roman [TT] found the result for null projection. Where we 
have 1/24, they had {v ■ £)^/12, where v is the unit tangent vector to the path of integration. 
Here v ■ i = = l/(5\/2), from Eq. fl5.44p . so the results agree. 

The remaining term requires more attention, because of the r dependence in yi. We 
write 

( 5 . 114 ) 


with 


Then 


Ri = 

Gi(r) = f 


— I drGi(r)^e 


) tt _ _ 
di{yiit,T))g (t- 

) 

= ^Gr(O) 


9[t + 


(5.115) 

(5.116) 


Applying the r derivatives to Gi gives 


GTir) 


dt 


T = 0 




r=0 


^(^^ + 3^2/1 (t,r) 




T=0 


- 4g"'(f)9(l) + 3g''(t) ) 


( 5 . 117 ) 


where the terms with an odd number of derivatives of the product of the sampling functions 
vanish after taking r = 0. 

Now yi depends on r and i only through t” = t + {\ — l/2)r, so using Eq. (I5.67p . we can 
write 

^g.(f.r) = T /‘dAF,);") = . (5.118) 

Then we integrate by parts and put all the derivatives on the sampling functions g, 


B 


1 — 



2 r dX (x - H%'i(t)(3g"(i)" + 4g'(i)g"'(i) + 9{i)g'"'(t)) 
+3 dX (a - 0 ' F(*)(g""(i)g(i) - g"(i)t 
+lym9""mi) - 4g"'(()g'(*) + 3g'W) ■ 


(5.119) 


Since we set r = 0, Yi has no A dependence and we can perform the integral. The result is 

1 

= dm{t){g"{t)^ - 2g"\t)g'{f) + 2g""{^g {^). (5.120) 
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For the term proportional to r_^, we have 


OO J/- /•CO 


B2= I drG2(r)4e“'^"- 

-/n "7F ./_ ,-v^ T*_ 


(5.121) 


where 


G^ir) = 


dty 2 {t,r)g g + . 


We calculate this Fourier transform in Appendix IB.21 and the result is 

1 


Applying the derivatives to G 2 gives 


G'"(r) 


dt 


T = 0 -OO 




r=0 


B 2 = gG'2"(0). 






T=0 


(5.122) 


(5.123) 


(5.124) 


Again the only dependence of 1/2 on r is in the form of t" so we can integrate by parts 


B. = 


noo rl 

/ dt d\ 
1-00 Jo 


2 ( A - - ) Y2m3g'{t)g"{t) + g{t)g"\t)) 


+F U - F y 2 {B){g''\t)g{t) - g''{t)g\t)) 


and perform the A integrals 


(5.125) 


B 2 = — I dtY 2 {i){g'{t)g'\t) - 3g''\Yjg{^). 


— OO 
-2 


For the term proportional to r_ , we have 


(5.126) 


OO /•OO 


B3= I - I drG'3(r)4e-'«A 

./n vr /__ T1 


(5.127) 


where 


GYt) = 


dt{ci{i) + ysit, T))g 9 (t + 


The result from Sec. 13.41 is 


Bs = ^G'm 


Applying the derivatives to G 3 gives 


G'^{t) 


dt 


T = 0 J -OO 


(GJ , 


9^^ + + y 3 {t}){g'\i)g{i) - g'ii)"^^ 


T=0 


As before, we integrate by parts 


(5.128) 

(5.129) 

(5.130) 


^00 /•! 


^3 = X / dt dX 

^ J-00 Jo 


2 (^A - -J (1 - X)XY3{t}{g'{i)^ + 9{^9'\i)) 

+l(fi(«) + (1 - A)Ay3(r))(9"(()9(«) - g'(ff) 


(5.131) 
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Integrating in A gives 


Bs 





+ ^y3{^{^g''{t)g{t) 


2g'm 


The three remaining terms have Fourier transforms given in Sec. 13.41 


B 4 

B, 

Be 


where we added 



dTg'(t + ^)9{t-^) 


dtg{i)‘^{c3{i) + C5{i)) 


In |r|c 2 (t) sgnr 




(5.132) 

(5.133a) 

(5.133b) 

(5.133c) 


C5(^ = -(2a + b)Riiu{i ), (5.134) 

which is the local curvature term from Eq. (I5.42p . 

The bound is now given by Eqs. fl5TT^ . fl5T13ll . flTmill . flTM . fl5T^ . fl5T^ . 


5.6 The inequality 

We would like to bound the correction terms Bi through Bq using bounds on the curvature 
and its derivatives. Using Eq. fl5.2p in Eq. fl5.68ap . we hnd 

|Fi(f)| < • (5.135) 

27r^ 

We can use Eq. fl5.135p in Eq. 05.1201) to get a bound on |i?i|. But will not be interested in 
specihc numerical factors, only the form of the quantities that appear in our bounds. So we 
will write 

l^il < (5.136) 

where [g] is an integral of some combination of the sampling function and its derivatives 
appearing in Eq. O5.120p . We will need many similar functionals Jn\g], which are listed at 
the end of the section. The number in the parenthesis shows the dimension of the integral, 

( 5 - 137 ) 

Similar analyses apply to B 2 and B^ and the results are 

1^21 < 4^^[5']|^“I-Rmax (5.138a) 

l^sl < 4^^[5']^max- (5.138b) 

Among the rest of the terms in B there are some components of the form Rated,u which 
diverge after boosting to the null geodesic, as shown in Ref. [23] • However we can show that 
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these derivatives are not a problem since we can integrate them by parts. Suppose we have a 
term of the form 

poo poo 

Bn= di dTLn{T,t)Rabcd,u{x) , (5.139) 

J OO j —OO 

where Ln{T,t) is a function that contains the sampling function g and its derivatives. The u 
derivative on the Riemann tensor can be written 


Bated,u B, 


abed.t 


-B. 


abed.v 


(5.140) 


The term can be reorganized the following way by grouping the terms with t and v,x,y 
derivatives 

poo poo 

Bn= dt dTLn{T,t){A'^^'^Rated,t{x) + Rated,a{x)) , (5.141) 

J OO J —OO 

where A'^'^'^'" are arrays with constant components and the subscript n denotes the term they 
come from. Here the greek indices a,^, - ■ ■ = v,x,y. The term with one derivative on a can 
be bounded while the term with one derivative on t can be integrated by parts, 

poo poo 

Bn = - dt driL^i, r)Hf "R,i„,(x) + r)H“R,i„,,„(x)). (5.142) 

J OO J —OO 

where the primes denote derivatives with respect to t. The sampling function is so L'{t, t) 
is still smooth and the boundary terms vanish. Now it is possible to bound this term, 

poo pOO 

|B„|</ dt dr(|L:,(f,r)|ofB„„ + |L„(f,T)|oWB:„„), (5.143) 

J OO J —OO 


where we defined 


(m) 




(5.144) 


The same method can be applied with more than one u derivative. 

Now we apply this method to the integrals R 4 , B^ and Bq of Eq. fl5.133p . We start with 
Bi, which has the form 


R 4 = 


dr In |rI sgn r / dtL^lt, r) Rated,tt(x) + Rated,ott{x) 


. ( 5 . 145 ) 


Liit, t) = g{t + T/2)g'{t -rjl.). 


where 

After integration by parts 

poo poo / 

R4 = / drln|r|sgnr / di(L'l(t,T)Af‘'‘^Ratcd{x)-L[(t,T)Af‘'‘^°'Ratcd,t 


(5.146) 


[X 


+L4,{t,T)Af‘"^°'^Rated,ag{x)^ ■ 


(5.147) 
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Taking the bound gives 


IB4I < fl ^‘^"'isi-Rl.r.’x 


(5.148) 


Reorganizing based on the number of t derivatives gives 


B. = 


dtL,{t) I Rated, Rated,ai{x) + (x) ) (5.149) 


di{ L''{i)Af^^Rated{x) - L',{t)Af^^Rated,a{x) + L^{i)Rated,afd{x) ) , 


where 


and the bound is 


Finally the remainder term is 


^00 poo 


LS) = 


\B,\<^jt^\g]R^:^^. 


(5.150) 


(5.151) 


Re = dt drUilr) dn / d\( Af‘^'^i\n)Rated,tttiX^) 


' —00 J —00 


+yl“(A, n)Rabed,att{XX}) + kl“^(A, n)R 

abed^a /3t{Xfl) 

+^f'"'^"/^^(A, Q.)Rated,al3^{X^) | SgU T , 


(5.152) 


where we changed variables to A = r/r and now arrays have components that depend 

on A and f2, and 

LQ{i,T) = g{i-Tl 2 )g{i+T/ 2 ) . (5.153) 

After integration by parts 


POO poo 


' —OO O —OO 


B, = I dr I di I dn I dxn'^iT,t}Af‘^\\,n)Ratedixxi) 


+Le{T, t)M“(A, n)Rated,a{Xn) + L,{t, f)'A“^(A, n)Rat,e.g{Xn) 

+Lq{t, f)Af'''^"^^(A, il)Rat,a0'YiXXl) | SgU T . (5.154) 


We dehne constants 


[ dXA, 

.j ^ o J i/0 


abed ajd... 
6 


(A,fl) , 


(5.155) 


abed a0... 


and now we can take the bound 


IBel < E 4‘"”’[9lfl££. 


(5.156) 
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Putting everything together gives 


6 3 


(m) 

^max 


(5.157) 


' —OO 

poo 


’ —OO 

poo 


’ —OO 

poo 


’ —OO 

poo 


’ —OO 

poo 


' —OO 

poo 


dt 


’ —OO 

poo 


’ —OO 

poo 


’ —OO 

poo 


’ —OO 

poo 


’ —OO 

poo 


' —OO 

poo 


' —OO 

poo 


dt 

dt 


' —OO 

poo 


— OO 

poo 


' —OO j —OO 


(5.158a) 

(5.158b) 

(5.158c) 


5 < r 2 ( i?o + 4^-") [g ]+ E E 

^ n=l n=4 m=0 

The functionals Jn\g] are 

OO 

dt{a^,\g""{t)\g{t) + ai 2 |/'(t)^'(t)| + a,,g"{tf 

) 

dt{a2i\g”\t)\g{t) + a22\g'\t)g\t)\) 
dt{a3i\g''{t)\g{t) + a 32 g'{tf) 

dt I dt'\\Yi\t-t'\\{a4,i\g'"(t')\g(t) + a42\g"(t)g'(t')\) (5.158d) 

) 

dt I dt'\ln\t - t'\\{a43\g''{t')\g{t) + aM\g'{t)g'{t')\) (5.158e) 

) 

dt’ |ln \t - t'\\a 45 \g'{t')\g{t) 

) 

dt{a5i\g"{t)\g{t) + a 32 g'{tf) 

) 

dt a53\g'{t)\g{t) 
dt a^^gitY 

dt I dt'Mg''\tMt') + a, 2 \g'\t)g'{t')\) 

) 

dt I dt'{aG3\g''{t)\g{t') + aG4\g'it)g'{t')\) 

dt’ a65\g’{t)\g{t’) 
dt’ amg{t)g{t '), 


= 

4"’la] = 
4‘’[sl = 
4‘’[9l = 
4"’[9l = 
4'‘’[9l = 

4‘’[9l = 
4“’[9l = 

Jt'la] = 

4‘’[9l = 
4“’l9l = 
4'‘’|9] = 

4^''’|9] = 


(5.158f) 

(5.158g) 

(5.158h) 

(5.1581) 
(5.158j) 
(5.158k) 

(5.1581) 
(5.158m) 


where ank are positive constants that may depend on 

We can change the argument of the sampling function, writing g{t) = f{t/to), where / is 
dehned in Sec. 15.21 and normalized according to Eq. (15.111) . so Eq. (I5.157p becomes 


dtTuu{w{t))g{tf > 


9-2 


^0 


r-to 


2A'kHi 


0 J —tQ 


dtf”{t/toY + E 


n=l 


6 3 


n=4 m=0 ' 


(5.159) 
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where we used Jn\g] = We can simplify the inequality by defining 



[ f'iafda = f'\t/tofdt , 

(5.160) 


6 

71=4 

(5.161) 

and 



(5.162) 

Then Eq. fl5.159p becomes 

j dtT^u{w{t))g{tf > 


(5.163) 


r-2 f 1 ^ ^ 1 

\ +E + E . 

0 L m=0 n=l ) 


We will use this result to prove the achronal ANEC. 


5.7 The proof of the theorem 


We use Eq. fl5.163p with w{t) = <hy(? 7 , f) and integrate in rj to get 




Tuu{^v{v,'t))f{t/tof > 


Vo 

6Hl 


247r2 




X 


u\3—n 


F(^)r 


max^O 


n—1 


m=0 


71=1 


(5.164) 


As (5 —)■ cxo, fo —t 0 but F‘'^\ F^'^\ Rmax, and R^2c are constant. Now = x'^/6, and using 
Eqs. fl5.10p . fl5.16p . |x“| < ui + \/26to. Thus as 5 —)■ oo, —)■ 0. Therefore only the hrst term 
in braces in Eq. fl5.164p survives, so the bound goes to zero as 

^ ~ , (5.165) 

Equation fl5.164p is a lower bound. It says that its left-hand side can be no more negative 
than the bound, which declines as But Eq. fl5.19|) gives an upper bound on the same 

quantity, saying that it must be more negative than —Ato/2, which goes to zero as R ~ 5““. 
Since a < 1/3, the lower bound goes to zero more rapidly, and therefore for sufficiently large 6, 
the lower bound will be closer to zero than the upper bound, and the two inequalities cannot 
be satished at the same time. This contradiction proves Theorem 1. 

The ambiguous local curvature terms do not contribute in the limit t]q ^ oo because they 
are total derivatives proportional to 

rm 

/ = 0 . 


-m 


(5.166) 





















Chapter 6 
Conclusions 


In this thesis we presented the derivation of quantum inequalities and a proof of the averaged 
null energy condition in curved spacetimes. Using a general quantum inequality derived by 
Fewster and Smith na, we first derived a bound for a quantum inequality in flat spacetime 
with a background potential, a case with similarities to the curved spacetime case, 


[ dtg{tf{T[r)Ut,0)> 

Jr 



h + -UnaxJ2 + U: 


11 


+u 


w 

max 


ll 7 r + 1 
IGtt 


Jr + 


' 2 '^^ J ( o7 J 1^1 

27r + 1 




24 

(dJg + J 7 ) 


Ja 

( 6 . 1 ) 


where the Jj integrals are given in Eq. fl3.117p . We then calculated the bound for a timelike 
projected quantum inequality in curved spacetime. 




IGtt^ 


U 4“ „ Rraax. J 2 

6 


( 6 . 2 ) 


+ R. 


f^ + 167r^(24|a| + 11|&|) ) Ja 


R 


n 

““ 60^" ' V40 
1637r + 14 


/// 


967r 




where the Jj’s in these case are given in Eq. fl4.106p . The importance and application of these 
results, for example in the special case of vacuum spacetimes, is discussed in Chapters [3] and 01 
Next we presented the derivation of a null projected quantum inequality in curved spacetime 


^ n=l n=4 m=0 

(6.3) 

where integrals are presented in Eq. (I5.158p . Finally we used this result to prove achronal 
ANEC in spacetimes with curvature. 



TabttdX > 0 . (6.4) 

As discussed in the introduction. Ref. [ 21 ] showed that to have an exotic spacetime there 
would have to be violation of ANEC on achronal geodesics, generated by a state of quantum 
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fields in that same spacetime. The result discussed above concerns integrals of the stress-energy 
tensor of a quantum held in a background spacetime; we have so far not been concerned about 
the back-reaction of the stress-energy tensor on the spacetime curvature. Thus we have shown 
that no spacetime that obeys NEC can be perturbed by a minimally-coupled quantum scalar 
held into one which violates achronal ANEC. This analysis is correct in the case where the 
quantum held under consideration produces only a small perturbation of the spacetime. Thus 
no such perturbation of a classical spacetime would allow wormholes, superluminal travel, or 
construction of time machines. 

What possibilities remain for the generation of such exotic phenomena? 

Could it be that a single ehect both violates ANEC and produces the curvature that 
allows ANEC to be violated? The following heuristic argument casts doubt on this possibility. 
Suppose ANEC violation and NEC violation have the same source. We will say that they are 
produced by an exotic stress-energy tensor Texotic- This Texotic gives rise to an exotic Einstein 
curvature tensor, 

exotic ^Planck ^exotic • (6.5) 


It is Gexotic that permits Texotic to arise from the quantum held. Without Gexotic, the spacetime 
would obey the null convergence condition, and so, since Texotic violates ANEC, it would have 
to vanish. A reasonable conjecture is that as Gexotic —^ 0, Texotic —)■ 0 at least linearlyji] We can 
then write schematically 

|Texotic| ^ I ^l^exoticl , (6.6) 


where I is a constant length obeying I ^ /pianck- The parameter /, needed on dimensional 
grounds, might be the wavelength of some excited modes of the quantum held. Equation fl6.6p 
is schematic because we have not said anything about the places at which these tensors should 
be compared, or in what coordinate system they should be measured. Combining Eqs. (I6.5p 
and fl6.6p . we hnd 


I .^exotic I ^ (^Planck 


/I) 


^IT -I 

I exotic I 


(6.7) 


which is impossible since I S> /pianck- 

Given the assumptions of this work, it appears that the only remaining possibility for 
self-consistent achronal ANEC violation using minimally coupled free helds is to have hrst a 
quantum held that violates NEC but obeys ANEC, and then a second quantum held (or a 
second, weaker ehect produced by the same held) that violates ANEC when propagating in 
the spacetime generated by the hrst held. The stress-energy tensor of the second held would 
be a small correction to that of the hrst, but this correction might lead to ANEC violation on 
geodesics that were achronal (and thus obeyed ANEC only marginally), taking into account 
only the hrst held. 

There is also the possibility of diherent helds. We have not studied higher-spin helds, but 
these typically obey the same energy conditions as minimally-coupled scalars. 

If one considers quantum scalar helds with non-minimal curvature coupling, the situation 
is rather diherent. Even classical non-minimally coupled scalar helds can violate ANEC inEi, 
with large enough (Planck-scale) held values. However, as the held values increase toward 
such levels, the ehective Newton’s constant hrst diverges and then becomes negative. Recently 
an even stronger result has been proven; the ehective Newton’s constant has to change sign 
between the two asymptotic regions on diherent ends of the wormhole [5] . Such situations may 
not be physically realizable. If one excludes such held values, some restrictions are known, 

^Not, for example, changing discontinuously for infinitesimal but nonzero Gexotic or going as 
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but there are no quantum inequalities of the usual sort |T2l US], and there are general |3l] 
and specihc [HU [32] cases where conformally coupled quantum scalar helds violate ANEC in 
curved space. It may be possible to control such situations by considering only cases where a 
spacetime is produced self-consistently by fields propagating in that spacetime, but the status 
of this “self-consistent achronal ANEC” for non-minimally coupled scalar fields outside the 
large-held region is not known. 
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Appendix A 

Multi-step Fermi coordinates 


In this Appendix, we generalize the Fermi coordinate construction to allow an arbitrary number 
of arbitrary subspaces (and an arbitrary number of dimensions d), rather than just a timelike 
geodesic and the perpendicular space. First, we construct the generalized coordinate system, 
then we compute the connection, and finally we compute the metric in the generalized Fermi 
coordinates. Results shown in this Appendix are used throughout the main body of the thesis. 
Here Greek indices a, jS ... refer to tetrad components while latin a,b,... to coordinate basis. 
This notation is not used in the rest of the thesis. 


A.l Multi-step Fermi coordinates 

Consider a d-dimensional Riemannian or Lorentzian manifold {J\A,g). We will start our con¬ 
struction by choosing a base point p E A4. We decompose the tangent space Tp into n 
subspaces, Tp = A^p'^ X A^p^ X Af ^... X A^p'^ so that any V E Tp can be uniquely written as 
V = V(i) -|- V( 2 ) -1- V^s) V(„). We choose, as a basis for Tp, d linearly independent vectors 

{E(^a)} adapted to the decomposition of Tp so that for each m = 1... n, {F(q) |q! E Cm} is a basis 
for where Ci, C 2 ,... c„ is an ordered partition of {1... d}. Thus each V(m) = 

The vectors {F(q,)} need not be normalized or orthogonal. 

The point corresponding to coordinates x°' is then found by starting from p and going 
along the geodesic whose whose tangent vector is V(i), parallel transporting the rest of the 
vectors, then along the geodesic whose tangent vector is V( 2 ), and so on. An example is shown 
in Fig. lA.ll With that general construction of multi-step Fermi coordinates we can define a 



Figure A.l: Construction of 3-step Fermi coordinates in 3 dimensions. We travel hrst in the 
direction of F'(i), then F( 2 ), then F( 3 ), parallel transporting the triad as we go. 
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general Fermi mapping q = Fermip(l/) given by 


m = p 

g(i) = expp(V(i)) 

g(2) = expg^^j(l/(2)) (A.l) 

q = q{n) = exp^^_^j (!/(„)) 

From that general constrnction we can retnrn to the original Fermi case by choosing ci = {t} 
and C 2 = {x,y,z}. In the Lorentzian case, we conld also choose a psendo-orthonormal tetrad 
Eu, Ex, Ey, with Eu and Ey nnll, Ey ■ Ey = — 1 , and other inner prodncts vanishing, and 
Co = {m} and ci = {v,x,y}. 

For later use we will define 

m 

V(<„, = Y. = E(^-2) 

«GclUC2U'"UCm ( = 1 

Then we can write q^m) = Fermip(V(<m))- 


A.2 Connection 


We will parallel transport our orthonormal basis vectors F^(q,) along the geodesics that generate 
the coordinates, and use them as a basis for vectors and tensors throughout the region of Ai 
covered by our coordinates. Components in this basis will be denoted by Greek indices. We will 
use Latin letters from the beginning of the alphabet to denote indices in the Fermi coordinate 
basis. Of course at p, there is no difference between these bases. 

Latin letters from the middle of the alphabet will denote the subspaces of Tp or equivalently 
the steps of the Fermi mapping process. 

We would like to calculate the covariant derivatives of the basis vectors, V, which are 
connected with connection one-forms, see for example Ref. [39] by 


^pa5 V'yS^P^(a) ’ 


(A.4) 


because we are using a orthonormal tetrad basis. 

We can then calculate the covariant derivative of any vector field V 
/(A) as 


DV^ 

dX 


dVf^ 

~dX 





V„E: 


(y) 


V^Eyj along a curve 


(A.5) 


To evaluate VpEy) at some point qi = expp(A), consider an infinitesimally separated 
point q 2 = expp(A -|- E(^y)dx). The covariant derivative of Ey) at qi is the difference between 
E(a){(l 2 ) parallel transported to qi and the actual Ey'^^qi), divided by dx. That difference is 
the same as the change in Ey^ by parallel transport around a loop following the geodesics 
from qi backward to p, the infinitesimally different geodesics forward from p to q 2 , and the 
infinitesimal distance back to qi. We can write this loop parallel transport as an integral over 
the Riemann tensor. 
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^2 



Figure A.2: The covariant derivative 'SJpEi^a) is the change in E(^a) under parallel transport 
along the path qi ^ p ^ q 2 ^ Qi, divided by dx. The parallel transport can be decomposed 
into a series of transports clockwise around trapezoidal regions with sides XE^p^dx and XdX. 


Let us hrst consider the Riemannian case, as shown in Fig. IA.2l The total parallel transport 
can be written as the sum of parallel transport around a succession of small trapezoidal regions 
whose sides are XE{_p) and XdX. By using the definition of the Riemann tensor we have 

XpEj. = f dX R\sp{XX)XX ^. (A.6) 

Jo 

Here R is evaluated at the point expp(AX), which we have denoted merely AX for compactness. 

Equation flA.6j) reproduces Eq. (13) of Ref. [29]. Note, however, that Eq. flA.6|) is exact 
and does not require R to be smooth, whereas that of Ref. |29| was given as first order in R 
and was derived by means of a Taylor series. 

We see immediately that the covariant derivative of any E(^a) at iu the direction of X 
vanishes. This happens simply because changes with dX in the direction of X correspond to 
additional parallel transport of E(q,). 

Let us now consider the general case where there are n steps, and compute XpE(^a)- Since 
the coordinates are adapted to our construction, the index (3 must be in some specific set 
Cm, which is to say that the direction of the covariant derivative, E(^), is part of step m in 
the Fermi coordinate process. We will write the function that gives that m as m{/3). Some 
particular cases are shown in Fig. IA.3I 

If m = n (leftmost in Fig. IA.3|) . only the last step is modified. The integration is exactly 
as shown in Fig. IA.21 except that it covers only the hnal geodesic from X(<„) to X(„), 

XpEj.= [ dXR\sp{X^^n) + XX^n))XXf^y (A.7) 

Jo 

If m < n, then we are modifying some intermediate step, and the path followed at later 
steps is displaced parallel to itself. In that case we get an integral over rectangular rather than 
trapezoidal regions, as shown in Fig. IA.41 For general m there is a contribution for each step 
j > m. The j = m contribution integrates over trapezoids that grow with A, while the j > m 
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Figure A.3: Original and displaced geodesics for Fermi coordinates with n = 3 and m = 3, 2, 
and 1. 


^2 



E(p)dx 


Figure A.4: Part of the calculation of V/ 3 -F(a) in the case n = 2, m{/3) = 1. The geodesic of 
the first step has been modified, causing the geodesic in the second step to be displaced. The 
parallel transport integrates the Riemann tensor over a series of rectangular regions between 
the 2 second-step geodesics. 
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contributions integrate over rectangles with fixed width dx. We can write the complete result 




+ XX(^,))Xl 


J) 


j=m 


where m = m{/3) and 


^jm 


(A) = 


1 j ^ rn 
X j = m . 


(A.8) 


(A.9) 


Equation flA.Sj) is exact and includes Eqs. flA.6IIA.7p as special cases. 

Consider the case where Ci consists only of one index. If m > 1, there is no j = 1 term in 
Eq. flA.Sp . If m = 1, then fd is the single index in ci, and = 0 unless 6 = (d, so the j = I 
term in vanishes because Ra^yS/s is antisymmetric under S ^ /d. Thus there is never a j = 1 
contribution to Eq. (lA.Sp when there is only one index in Ci. 

Now suppose X lies on the first generating geodesic, so X(j) = 0 for j > 1. Then all j > 1 
terms vanish in Eq. (lA.Sp . So if ci consists only of one index, all Christoffel symbols vanish at 
X. This is well known in the case of the usual Fermi coordinates. 


A.3 Metric 


Now we would like to compute the metric g at some point X. Specifically, we would like to 
compute the metric component Qah in our generalized Fermi coordinates. 

We will start by considering the vectors Zf^a) = d/dx^. These are the basis vectors of the 
Fermi coordinate basis for the tangent space, so the metric is given by gab = Z(^a) • Thus 
if we compute the orthonormal basis components we can write gab = da^^{a)^{b)- 

Again we will start with the case of Riemann normal coordinates. Let W{t, s) be the point 
exppS{X+tE^a))- Define Y = dW/dt and V = dW/ds. Then Y{X) = Z^a) and 
The components of Z(^a) at X can be calculated by integration. 


ZUX) = Y^X) = 


'0 


, dY^{sX) 

ds ---. 

os 


Because the orthonormal basis is parallel transported we can write 


±Yl> = —. 

ds ds 

By construction, the Lie derivative LyY = 0 and thus [22l Ch. 4] 

DY _ DV 
ds dt 


From Eq. flA.Sj) we have 


(A.IO) 


(A.ll) 


(A.12) 


DV^ 

dt 


dV^ 

dt 




if + + O(R^). 


(A.13) 


where we have retained instead of writing to make it clear that the covariarrt 

derivative is with respect to the orthonormal basis. 
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From Eq. flA.6jl we have 


V„Ej^^)(sX)= / dXXR^^sa{XsX)sX^ = - I dX XR^^5a{XX)X^ 


(A.14) 


Taking t = 0, E is just X. Combining Eqs. flA.10IIA.14|) . we find 


4,(A-) = 


ds 


5f + (5“ / dAAi?%„(AX)X^X' 


+ 0{R^ 


= 5f + (5“ / dXX{l - X)R^^5a{XX)X^X^ + 0{R^). 


(A.15) 


From Eq. flA.15ll . the metric is given by 


gab = Vab + 26:6^ta [ dX A(1 - X)R^^sp{XX)X^X^ + 0{R^). (A.16) 

Jo 

Equation flA.lbp reproduces Eq. (14) of Ref. 

Next let us consider the case where there are n steps in our procedure. We will define a 
set of functions Wj as 

Wj{s) = Fermip(X(<j) + sX^j)). (A.17) 

The path Wj{s),j = 1.. .n,s = 0...1 traces the geodesics generating the Fermi coordinates 
for the point X. Now consider Z(^a) = djdx'^. Let m = m{a), so Z(^a){p) G Then let 


{ sX(j) j <m 

X«j) + s(X(j) + tE(^a)) j = m 
X«j) + tE(^a) + sX(j) j >m 


(A.18) 


Let Y = dW/dt and V = dW/ds as before. To find Zi^a) we now must integrate over a 
multi-step path from p, 

aY>{wM) 


4)W = E /'* 

j=l -'0 


ds 


The generalized version of Eq. flA.lSp is 


(A.19) 


0 J ^ 

DV^{W,{s,t)) ^ ^ =\s! + + 0{R^) j = m (A.20) 


dt 


dt 


(7) 

+ 0(R2) 


j > m. 


Now 



dX Qjj^jYiiX^R 'ySa(^X(^^if;'j Y XX^i^'j'jX^f^'j 


(A.21) 


^Ref. [29] uses the same sign convention for as the present thesis, but the opposite convention for gab 

and consequently also for Rap-yS- 
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where 


= i (A-22) 

ys k = j. 

The k = j term is analogous to Eq. flA.14D . while the others have no dependence on s. 
Combining Eqs. flA.llllA.12IIA.1911A.2ip we get 

Zf.)(A') = if + q' + 0(ii^) (A.23) 

where 

^ ri rskj{s) 

L ds / d\akm{^)R^ 

j=m k=m ^ 

n j „i 

d\ akm{\)hkM)R^-,UXi<k) + (A.24) 

j=m k=m 


where m = m{a) and 

h (W 

(A.25) 

Thus the metric is 

Qab = VagZ^a)F^b) ~ dab + Fab + Fba + 0{R^) 

(A.26) 

where 

Tl 1 1 


Fag = 

'fc J pi 

j=mk=m ^ 

(A.27) 

where m = m{/3). 




Florides and Synge [16] construct coordinates by taking geodesics perpendicular to an 
embedded submanifold. Our construction for n = 2 is of this kind in the case where all basis 
vectors with indices in ci he tangent to the surface generated by all first-step geodesics. This 
will be so if Ra-ysp = 0 everywhere on this surface whenever ^,6,(3 G ci and a G C 2 . In that 
case, Eq. flA.27D agrees with Theorem I of Ref. |16j . 

Now, consider again the case where Ci contains only one index. As discussed with respect 
to Eq. flA.Sp . ii (3 G Ci, there is no nonvanishing k = 1 term in Eq. flA.27D . Thus gab = Vab 
everywhere on the hrst generating geodesic. This is also well known in the usual Fermi case. 

Now suppose Cl consists only of one index and furthermore n = 2. The only possible term 
in Eq. flA.27H is then j = k = 2, so 


Fag — / dAa2m(A)(l — A)RQ,.y5/3(X(i)-|-AX(2))X(2)X^"^2)- 


(A.28) 


where m = m{(3). Equation Eq. flA.28p is equivalent to Eq. (28) in Ref. 
the generating curve of the Fermi coordinates is a geodesic. 


in the case where 
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A.4 Regularity of the coordinates 

Riemann normal coordinates cannot in general be defined over the entirety of a manifold, 
because there might be points conjugate to the base point p. At such a point, some infinitesimal 
change to the coordinates would yield no change to the resulting point, and the metric would 
be singular. 

Similar considerations apply to Fermi normal coordinates [27]. No trouble can occur in the 
first step, because that consists merely of traveling down a geodesic. Along the geodesic, and 
thus by continuity in some neighborhood surrounding the geodesic, Fermi normal coordinates 
are regular. If we attempt to extend beyond this neighborhood, we may find points that are 
conjugate to the generating geodesic. In such places the metric will become singular. 

The situation here is more complicated. The metric will be singular whenever an infinites¬ 
imal change in coordinates fails to yield a change in the location of the resulting point. But 
when there are more than two steps, the result can no longer be described in terms of conjugate 
points. Nevertheless, it is easy to see that one if one chooses a sufficiently small neighborhood 
around p, all multi-step Fermi coordinates will be well defined, since any such coordinates 
approach Riemann normal coordinates when all coordinate values are sufficiently small. 

In case Ci contains only one index, multi-step Fermi coordinates will be well defined in a 
neighborhood of the initial geodesic, because when all coordinate values except for the first 
are small, the multi-step Fermi coordinates approach the regular Fermi coordinates. There is 
no particular advantage to having a single index in any later Cm- 

One can get a simple condition sufficient for the existence of multi-step Fermi coordinates 
in a small region by looking at Eqs. flA.26p and flA.27p . As long as \Fab\ <C 1, the metric Pab 
cannot degenerate. Thus the coordinates will be well defined if [29] 

\Ra^f}s\{X^f I (A.29) 

throughout the region of interest, for all a, 7 , f3, 6, e. In the case where there is only one index 
in Cl, there is no contribution to Fab from X(i). Then it is sufficient for Eq. flA.29p to hold for 
e > 1. Thus if the first step is one-dimensional, it can be arbitrarily long [27], as discussed 
above. 
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Appendix B 

Fourier transforms of some 
distributions 

B.l Fourier transforms of some distributions involving 
logarithms 


In this appendix will compute the Fourier transforms of the distributions given by 

u{t) = ln|r| 
v{t) = ln(—r^). 

We write m as a distributional limit, 


where 


so its Fourier transform is 


u. 


u = lim Ue , 

e^0+ 


:(t) = In |r|e , 


Ue{k)= f drln|r|e-^l"le*^" = 2Re / drlnr =-2Re ^ ^ 


e — ik 


Thus the action of h on a test function / is 


u[f] = -2 lim Re [ ^^ 7 + ln(e 

e-i^o+ e — ik 


The term involving 7 is 


27 lim 

e^0+ 


dk 


k"^ + e 


:f{k) = -27r7/(0). 


In the other term we integrate by parts, 

ln(e — ik) 


— 2 lim Re / dk- 

£-> 0 + / „ 


(B.l) 

(B.2) 

(B.3) 

(B.4) 

(B.5) 

(B.6) 

(B.7) 


f{k) = — lim Im / dk f'{k)[ln{€ — ik)]‘^ 

e-ik e^o+ 

/ OO 

dk f'{k)[\n \k\ — i(7r/2) sgn A;]^ 

-OO 

/ OO 

dk f'{k) In \k\ sgn/c, 


(B. 8 ) 
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and thus 


■hi/] = tt / dfc/'(/c) In |/c| sgn/c — 27r7/(0). 

J —OO 

Since the Fourier transform of the constant 7 is just 27r'j6{k), the transform of 

w{t) = In |r| + 7 

is just 


■^[Z] = ^ / dfc/'(/c) In |/c| sgn/c. 


v(t) = limln(—(r — ieV) = 2 In |r| + 7ii sgnr . 


Now 


The Fourier transform of sgn acts on / as 

2iP f dk^ = -2if dkf{k)\n\k\, 

J —00 ^ J —00 

Putting Eqs. flB.9IIB.13p in Eq. flB.12p gives 

poo 

v[f] = Air / dk f'{k)\n\k\ — An'yf{0). 


(B.9) 

(B.IO) 

(B.ll) 

(B.12) 

(B.13) 

(B.14) 


B.2 Fourier transform of distribution r 


We follow the procedure of Sec. 13.41 to calculate 

B2= [ —[ drG2(r)s2(r)e-*«T 

Jo J—00 


(B.15) 


where 

G2(r) = j diy2{t, T)g (Z- ^ (Z+ • (B.16) 

and 

= (B.17) 

This is the Fourier transform of a product so we can write it as a convolution. The function 
S 2 is real and odd, so its Fourier transform is imaginary, but G 2 is also real and odd, thus the 
Fourier transform of their product is real. We have 


1 poo poo 

^^ = ^1 dCG2{-^ - OUC) ■ (B.18) 

We can change the order of integrals and change variables to rj = — ( which gives 

1 poo poo 

B 2 = dcj dvG2iv)s2{C) 

1 r°° ^ n 

= J dr]G2{v) J dCs2{C). (B.19) 
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The Fourier transform of S 2 is [20] 



S2(C) = 

-i7rC^0(C) , 

(B.20) 

and 


^ dC{-i7rC) = -y77^0(r/). 

(B. 21 ) 

From Eq. (IB. 19) we have 






1 dr] G 2 {r])r]^ . 

0 

(B.22) 

Using /'(O = we get 

j 

noo 

dvG'^iv). 

0 

(B.23) 


The function G 2 is odd but with three derivatives it becomes even, so we can extend the 
intergal 

^2 = y dh = gG"'(0). (B.24) 
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